e  School 


PARAMETER  PLANE  ANALYSIS  OF  SAMPLED  DATA  SYSTEMS 
(WITH  EXTENSIONS  FOR  CONTINUOUS  SYSTEM  APPLICABILITY) 

***** 
Aloysius  R.  Miller 


PARAMETER  PLANE  ANALYSIS  OF  SAMPLED  DATA  SYSTEMS 
(WITH  EXTENSIONS  FOR  CONTINUOUS  SYSTEM  APPLICABILITY) 


by 

Aloysius  R.  Miller 
ft 

Lieutenant,  United  States  Navy 


Submitted  in  partial  fulfillment  of 
the  requirements  for  the  degree  of 

DOCTOR  OF  PHILOSOPHY 

IN 

ELECTRICAL  ENGINEERING 

United  States  Naval  Postgraduate  School 
Monterey,  California 

19  6  5 


- 


PARAMETER  PLANE  ANALYSIS  OF  SAMPLED  DATA  SYSTEMS 
(WITH  EXTENSIONS  FOR  CONTINUOUS  SYSTEM  APPLICABILITY) 

by 
Aloysius  R.  Miller 

This  work  is  accepted  as  fulfilling 

the  dissertation  requirements  for  the  degree  of 

DOCTOR  OF  PHILOSOPHY 

IN 

ELECTRICAL  ENGINEERING 

from  the 

United  States  Naval  Postgraduate  School 


ABSTRACT 
Many  of  the  advances  in  servo-system  design  fostered  by  the  general 
availability  of  digital  computers  are  not  utilized  by  practicing  engineers 
because  they  discard  rather  than  complement  existing  methods.   An  exception 
is  the  parameter  plane  method  which  retains  the  principle  of  design  by  such 
familiar  specifications  as  damping  and  settling  time,  but  where  classical 
methods  display  the  effect  of  only  one  variable  on  these  specifications  the 
parameter  plane  method  can  display  the  simultaneous  effect  of  two.   In  this 
paper  the  restrictions  on  problems  suitable  for  parameter  plane  analysis 
are  eased  and  the  number  of  specifications  which  can  be  considered  on  the 
parameter  plane  is  enlarged  to  include  bltndwidth  and  steady  state  error. 
Parameter  plane  methods  are  also  adapted  to  dominant  mode  design  in  a 
manner  which  not  only  allows  selection  of  the  dominant  mode  but  simultan- 
eously guarantees  its  dominance.   These  Extensions  of  parameter  plane  theory 
are  related  to  both  continuous  and  sampled  data  systems. 
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CHAPTER  I 

In  recent  years  great  emphasis  has  been  placed  on  the  "controls  field" 
in  electrical  engineering.   With  this  emphasis  has  come  a  proliferation  of 
techniques  for  analysis  and  design  of  servomechanisms,  many  of  which  have 
not  yet  been  fully  exploited.   Among  such  techniques  are  the  parameter  plane 
methods  which  can  be  extended  in  both  applicability  and  in  specifications 
considered. 

Before  analyzing  the  parameter  plane  as  a  design  tool  it  may  be  wise  to 
establish  its  place  among  tools  presently  available  to  the  engineer.   A 
broadbrush  review  of  available  methods  should  be  sufficient  for  this  pur- 
pose.  Present  day  servomechanism  design  subdivides  conveniently  into  two 
disciplines,  direct  time  domain  analysis  and  frequency  or  transform  analysis. 
Of  the  two  the  former  appears  to  be  the  more  powerful  albeit  the  more  com- 
plex.  Design  by  this  procedure  ordinarily  includes  the  steps: 

a)  determination  of  signal  characteristics 

b)  selection  of  a  criterion  by  which  system  performance  can  be 
evaluated  (i.e.  integral  squared  error,  etc.) 

When  the  results  of  these  steps  are  clear  cut  and  definite  time  domain 
analysis  will  probably  be  the  best  tool  available  to  the  design  engineer. 
This  holds  particularly  for  sampled  data  systems  where  incorporation  of  a 
digital  computer  can  make  realizability  of  theoretically  obtained  results  a 
minor  problem.   Despite  its  power  time  domain  analysis  is  by  no  means  a 
panacea.   By  its  very  nature  the  time  domain  approach  demands  that  system 
structure,  its  inputs  and  its  objectives  be  well  defined.   These  are  by  no 
means  trivial  requirements  since  a  system  designed  on  the  basis  of  integral 


squared  error,  for  example,  need  not  have  acceptable  acceleration  character- 
istics, overshoot,  frequency  response,  etc.   Furthermore,  analysis  in  the 
time  domain  generally  provides  little  insight  into  possible  performance 
criteria  trade-offs  or  structural  changes  which  might  enhance  the  overall 
system. 

The  difficulties  of  obtaining  performance  characteristics  from  the 
integrodifferential  equations  which  describe  a  servo-system  have  to  a  large 
extent  been  overcome  by  the  use  of  transforms  which  reduce  the  operations 
of  differentiation  and  integration  to  operations  of  multiplication  and 
division.   When  transformed,  the  pole  zero  pattern  or  magnitude  phase  plot 
of  the  integrodifferential  equation  furnishes  an  indication  of  system  per- 
formance by  inspection.   In  addition  these  plots  often  indicate  system 
modifications  or  compensation  which  will  improve  system  performance.   The 
primary  disadvantage  of  transform  analysis  is  that  pole-zero  patterns, 
bandwidth  and  other  transform  criteria  abstract  the  time  domain  performance 
of  the  system  and  it  is  in  the  time  domain  that  the  system  must  ultimately 
perform  satisfactorily. 

A  great  number  of  servo  design  procedures  fall  within  the  category  of 
transform  techniques.   In  a  broad  sense  these  procedures  can  be  subdivided 
into  a  class  dominated  by  a  synthesis  concept  and  a  class  dominated  by  a 
compensation  concept  although  in  practice  neither  can  be  isolated  from  the 
other.   The  synthesis  approach  to  servo  design  is  accomplished  in  three 
steps. 

1)  determination  of  the  closed  loop  transfer  function  from 
specifications 

2)  determination  of  the  open  loop  transfer  function  from  the 
closed  loop  transfer  function  obtained  in  step  (1) 
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3)  synthesis  of  the  open  loop  transfer  function  with  due  regard 
to  components  which  are  not  at  the  engineer's  disposal. 
While  the  transformed  equations  are  sufficient  to  determine  the  system's 
time  domain  behavior  the  complex  relationship  between  the  two  thwarts  ready 
selection  of  the  "best"  transfer  function.   Hence  it  is  likely  that  the 
selected  transfer  function  which  meets  specifications  is  more  difficult  to 
synthesize  than  some  other  transfer  function  which  also  meets  specifications, 
Furthermore,  there  is  no  systematic  or  iterative  procedure  which  will  im- 
prove this  initially  selected  transfer  function.   To  this  liability  must 
be  added  the  problems  that  fixed  system  components  impose  on  the  choice  of 
possible  closed  loop  transfer  functions.   To  its  credit  the  synthesis  con- 
cept follows  a  logical  sequence  from  specification  to  final  realization. 
This  is  in  contrast  to  the  compensation  approach  which  attempts  to  "patch 
up"  an  unsatisfactory  initial  system  in  the  hope  it  can  eventually  be  made 
to  meet  specifications.   The  second  liability  has  also  been  partially  elimin- 
ated at  the  expense  of  simplicity.  With  the  parameters  at  his  disposal  the 
engineer  is  able  to  produce  the  system  which  most  closely  approximates  the 
one  initially  chosen.  [1]  However  this  closest  point  of  approach  is  not 
guaranteed  to  meet  specifications  nor  does  it  preclude  that  a  satisfactory 
system  could  be  obtained  with  available  parameters  were  a  different  choice 
for  the  "ideal"  closed  loop  transfer  function  made. 

Utilization  of  transform  techniques  to  design  compensation  for  an 
existing  skeleton  servo-system  is  an  approach  more  widely  employed  than  syn- 
thesis.  Two  steps  are  basic  to  design  by  compensation. 

1)  From  specifications  a  "best  choice"  for  the  fixed  elements  of 
the  servo-system  is  made 


2)  Then  compensation  for  these  elements  not  under  the  engineer's 
control  is  designed  to  insure  system  compliance  with  specifi- 
cations. 
These  steps  contrast  vividly  with  those  employed  in  synthesis,  largely  elim- 
inate the  drawbacks  of  that  method  and  offer  very  definite  advantages.   The 
realization  problem  is  eliminated  by  the  expedient  of  employing  only  realiz- 
able compensators  to  modify  the  fixed  skeleton  system.   Of  almost  equal  im- 
portance are  the  well  established  methods  available  for  design  of  compensa- 
tion.  They  are  for  the  most  part  simple,  powerful  and  provide  insight  in 
addition  to  numerical  solutions;  a  feature  which  is  almost  unique  to  these 
methods.   In  addition  systems  designed  from  a  compensation  concept  tend  to 
be  comparatively  inexpensive,  a  fortuitous  result  which  usually  occurs  when 
dealing  with  the  physical  compensation  rather  than  the  theoretical  system. 
However  design  by  compensation  also  has  its  disadvantages.   Most  serious  is 
ignorance  of  whether  the  possibility  for  compensation  exists;  the  fixed 
system  may  be  incapable  of  being  compensated  to  meet  specifications.   In 
this  respect  the  synthesis  concept  is  clearly  superior  for  whenever  an  accept- 
able theoretical  system  cannot  be  found  then  it  serves  no  purpose  to  attempt 
design  of  an  impossible  physical  system. 

The  preceding  thumbnail  sketches  suffice  to  acquaint  the  reader  with 
approaches  available  to  the  servo-system  designer.   However,  deficiencies 
inherent  to  an  approach  are  subject  to  the  "state  of  the  art"  and  to  the 
state  of  development  of  the  methods  employed.   The  compensation  approach 
illustrates  the  case  in  point.   In  an  abstract  formulation  its  merits  are 
not  outstanding  yet  the  simplicity  of  the  techniques  and  the  well  documented 
correlations  between  available  variables  and  design  criteria  serve  to  make 
this  the  most  prevalent  of  approaches.   Parameter  plane  techniques  also 
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advantageously  utilize  these  well  established  correlations  between  f requency 
response,  bandwidth,  pole-zero  patterns  and  time  domain  performance  specifica- 
tions such  as  rise  time,  maximum  overshoot  and  settling  time.   Furthermore 
the  essential  steady  state  nature  of  frequency  response  analysis  and  the 
one  parameter  limitation  of  root  locus  diagrams  are  largely  overcome  by 
parameter  plane  methods.   At  the  expense  of  simplicity  the  effect  of  simultan- 
eous variation  of  two  parameters  on  root  locations  can  be  displayed.   Band- 
width specifications  and  steady  state  error  requirements  can  also  be  indicat- 
ed on  this  plot.   Furthermore  such  troublesome  questions  as  second  order 
dominance  or  whether  dominance  is  possible  by  adjustment  of  two  parameters 
can  be  answered  by  inspection.   These  advantages  auger  more  widespread  applica- 
tion of  parameter  plane  methods  in  systems  of  up  to  moderate  complexity, 
especially  as  the  digital  computer  becomes  more  readily  available  to  compute 
the  required  curves. 

The  intent  of  this  chapter  was  to  introduce  the  parameter  plane  and 
establish  its  position  among  tools  available  to  the  engineer.   The  division 
of  servo-system  theory  into  various  approaches  was  arbitrary  although  it 
seemed  reasonable.   Similarly  critique  of  present  approaches  and  methods  was 
pursued  as  a  means  of  placing  the  parameter  plane  and  no  claim  for  complete- 
ness or  comprehensive  analysis  is  made. 


CHAPTER  II 

Concepts  of  the  Parameter  Plane 

The  characteristic  equation  which  is  obtained  by  equating  the  denominator 
of  the  closed  loop  transfer  function  to  zero,  is  probably  the  most  important 
single  factor  available  for  determining  system  performance.   In  deference  to 
convention  the  characteristic  equation  will  be  designated  by 

f(s)=  ^  aKS   =  °  (*> 

K=o 

where   s   is  the  complex  variable  of  the  LaPlace  transform 

n  is  the  order  of  the  differential  equation  describing  the  system 
a,  are  real  constants 
However  such  a  designation  is  an  over  simplification.   Consider  the  follow- 
ing points. 

1.  f(s)  which  is  an  equation  in  the  complex  domain,  can  properly 
be  regarded  as  two  simultaneous  real  equations.   Namely: 

Real  Part   [f  (S)|  =  ?R(S)  =  O 

Imaginary  Part   |_T  (S)J  =  tj  (•$)  =  O 

2.  The  complex  variable   s   is  in  reality  two  real  variables 
since  both  its  real  and  imaginary  parts  must  be  specified 
before   s   is  fixed.   Variables  commonly  used  to  specify  s 
are 

s  =  <r+  juj  =  u/N  [-  C  +  ij]/l-  tj 

For  parameter  plane  formulations  the  latter  specification 
is  the  more  convenient. 


3.   Invariably  a  servo  system  has  a  number  of  parameters  which 
can  readily  be  varied  to  suit  the  designer  (i.e.  system 
gain)  and  any  such  change  will  affect  the  coefficients  a.  . 
For  parameter  plane  purposes  it  is  convenient  to  assume  that 
the  £L  are  functions  of  two  such  variables  which  are  independ- 
ently at  the  designer's  discretion.  That  is:   a,  =  \(°(    »  /3  ) 
(The  existance  of  more  than  two  such  variables  does  not  in- 
validate parameter  plane  concepts  any  more  than  the  existance 
of  more  than  one  variable  invalidated  root  locus  concepts. 
Rather,  all  but  two  of  the  variable  parameters  should  be  fixed 
prior  to  any  single  application  of  parameter  plane  technique). 
With  these  points  in  mind  a  more  cogent  statement  of  the 
characteristic  equation  would  be 

k 


(  ©<  »  ft  f  £  »  UJM   )  -  Real 


2  &*(<*>  fi)S 


K=0 


fI  (  &  »  (3   »  £   »  H  )  *  Imaginary 


J   (2) 


2>i<*./D5' 


Equations  (2)  are  two  equations  in  four  unknowns.  Classically  the  coefficients, 
a,  ,  were  assumed  known  (the  designer  assigned  values  to  o(, 3   )  and  the  pro- 
blem was  reduced  to  one  of  finding  roots  for  the  characteristic  equation. 
If  these  roots  were  unsatisfactory  the  designer  would  assign  new  values  to 
o(   and yS  and  repeat  the  procedure.   This  procedure  is  somewhat  illogical  in 
that  unknowns  (variable  system  parameters)  are  initially  fixed  in  expecta- 
tion of  obtaining  knowns  (desired  root  locations).  A  better  procedure  would 
be  to  solve  for  unknowns  in  terms  of  the  knowns.   That  is,  solve  the  implicit 
equations  (2)  for  explicit  expression  of  0<  and  R   . 
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/a-iStf.uv)  <3> 

Direct  determination  of  system  parameters  in  terms  of  desired  s-plane  loca- 
tions £  and  U/N    is  now  feasible.   However,  s-plane  roots  are  not  generally 
restricted  to  a  single  point;  rather  requirements  are  such  as  to  place 
them  within  a  certain  region.   Suppose  that  a  number  of  s-plane  points  on 
the  boundary  of  this  region  are  chosen,  that  the  corresponding  values  of 
system  parameters  o<  and  j8    are  computed  using  equations  (3)  and  that  these 
values  are  plotted  as  points  in  a  cartesian  o(-j3   coordinate  system.  The 
curve  connecting  these  points  could  reasonably  be  expected  to  form  a  boundary 
for  the  o<-^0  plane  region  containing  values  of  system  parameters  which  place 
roots  of  the  characteristic  equation  within  the  desired  s-plane  region. 
(A  more  detailed  interpretation  of  c^-jS   plane  curves  is  reserved  for  the 
next  chapter).   Early  investigators  were  interested  in  this  concept  pri- 
marily as  an  aid  to  stability  analysis  and  therefore  restricted  their  atten- 
tions to  an  s-plane  contour  along  the  jw  axis  [4,5] .   The  subsequent  work  of 
Mitrovic  and  others  generalized  procedures  to  include  constant  £  contours 
but  were  still  limited  by  restricting  o(   and  ^8  to  be  the  last  two  co- 
efficients of  the  characteristic  equation  [3].  More  recently  Siljak  general- 
ized the  theory  to  allow  the  coefficients  of  the  characteristic  equation  to 
contain  any  linear  combination  of  o<,  and  G      [7],  This  paper  extends  the 
procedures  further  to  include  the  appearance  of  o<./Q    products  in  the  co- 
efficients.   Still  further  extension  to  appearance  of  o^c  and  /3   terms  and 
beyond  is  feasible  but  resulting  explicit  expressions  for  o(f£i<^)  and  /SC^^m) 
rapidly  become  prohibitive. 

Parameter  plane  (or  o(-S   )  curves  can  be  considered  from  an  alternative 
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viewpoint  which  is  more  suitable  to  mathematical  rigor.   Here  (  c<  ,  j3   )  and 
(  2  »  WH     )  are  two  sets  of  variables  which  are  related  by  either  implicit 
equations  (2)  or  explicit  equations  (3).   The  parameter  plane  (or  o(-j8 )   curves 
are  then  the  images  (or  maps)  of  the  chosen  s-plane  contours.   (The  reader  is 
cautioned  that  the  mapping  is  not  conformal  and  basic  relations  of  conformal 
mappings  are  not  preserved) . 

Parameter  Plane  Equations  (Linear  Coefficient  Case) 

This  formulation  is  an  abbreviated  version  of  that  advanced  by  Siljak 
and  as  such  presents  only  the  salient  points  [7].  Consider  the  character- 
istic equation  (1)  where  the  coefficients  a,  are  linear  functions  of  the 


parameters  o<  and  B 


,/s,s )  -  2  OxM-^  +  eJs*  = 


K=0 
In  this  equation  b,  ,  c,  and  e,  are  constants  and  s  is  the  complex  variable 

of  the  LaPlace  transform.   It  is  convenient  to  express  s  as 


*  Htf+a//-**]  (5) 


and 


s* 


=  u%  (j,  (- i)  + 1  y^F  wK  (-  e)|    k*o  (6) 


where  w  is  the  undamped  natural  frequency  and  £    is  the  damping  constant 
which  assumes  positive  values  in  the  second  quadrant  from  zero  on  the 
imaginary  axis  to  plus  one  on  the  negative  real  axis.  Observe  that  equations 
(5)  and  (6)  are  not  inconsistant;  rather  equation  (5)  serves  to  define  T,(-£) 
and  U,  (-  £  )  as  explicit  functions  of  £  . 

By  manipulation  of  equations  (5)  and  (6)  some  subsequently  useful  re- 
lations between  the  U,  's  and  T.  's  can  be  obtained. 

k       k 


Argument  relations 

Tk(-  £)  -  (-D*  Tk(5  )  (7) 

uk(-S  )  =  (-1)**'   uk(^  ) 

Recurrence  relations 

Tk+1(  5  )  "2  ^  Tk(  g  )  +  Tk-1(  £  )  «  0  (8) 

uk+i(  S>  ~2£  uk(£  )  +  ukl(  2)  =  0 

Interrelation 

V  *  }  =  ^  V  ?  >  "uk-i(  5  }  (9) 

Also  the  T,  (  £  )'s  an(*  U  (  ,£  )'s  are  found  to  be  Chebyshev  functions  of  the 
first  and  second  kind  respectively. 

Armed  with  these  relations  we  proceed  to  find  explicit  solutions  for 
o(  and  Q    .   Substitution  of  definition  (6)  into  the  characteristic  equa- 
tion as  expressed  in  (4)  yields: 

H 

k=o  l_        *  _l 

The  real  and  imaginary  parts  must  independently  equal  zero  to  satisfy  this 
equation. 


K*0 

Removing  common  factors  and  employing  relations  (7)  to  change  the  arguments 
of  the  U,  '  s  and  T  's  to  positive  zeta  results  in  the  expressions: 


K=0 


1"°  10 


Equations  (9)  which  express  the  interrelation  between  the  T,  's  and  U  's  can 

iv  K 


be  used  to  eliminate  T,  . 

k 


?  /Wm^ :aK(s) 


t-o 


]4«* 


«^fic*+t*)u£uje) 


=o 


g(-/)k  (d.  bK+  /3cK  +  e,)  utf  UK  ( £)]  =  o 


But  the  quantity  in  the  left  hand  bracket  is  identical  to  the  quantity  below 
it  which  is  equal  to  zero.   Incorporating  this  observation  allows  us  to  write 


2 

K=0 


(-/)    (o<bK+^c,-i-eK)H',U)(.1C5)  =  o 


(10) 


Equations  (10)  are  linear  in  o<  and  ft   and  are  easily  solved  for  their  explicit 
expression.   This  point  is  readily  apparent  when  equations  (10)  are  rewritten 
in  more  concise  notation. 


o(  B2(£,K,)  +■  /3Cg(?,uyv)  +  re({,ii4)»o 

where  by  definition 


(11) 


B2(  £  .  H 

cx(   £  .  ^ 

c2(  £  ,  ^ 
V  £  .  "* 

e2(  £  , h 


■1<- 

■2'- 


\)(t*s[)  "  -  ud 


)kLw„K  LL( s) 


D, 


>  C«Mf  "..<«)«! 


(12) 


)  cKw„*UK(S)   •     cx 


I    ■   -w 


- 
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In  the  interest  of  brevity  the  functional  dependence  of  B,  C  and  E  will  be 
omitted  when  such  omission  can  cause  no  confusion.  The  explicit  solutions 
f or  o<  and  3    obtained  from  equations  (11)  are: 

U  (  £  ,  <4  )  -   ^'5e~CJ^'  <13> 

B.Cg-BgC, 

/3  (  £  ,  14,  )  =   6a  g  -  6,  E> 

B,C2-  6BC, 

These  are  the  explicit  expressions  for  o^and  |3  which  allow  mappings  of 
points  or  contours  (except  for  points  on  the  real  axis)  from  the  s-plane 
into  the  parameter  (or  o(- B   )  plane.  Applicability  does  not  extend  to 
points  on  the  real  axis  because  the  derivation  included  division  by  a  factor 
of  (1  -  j»  )  which  is  identically  zero  thereon. 

Relevant  parameter  plane  equations  for  points  on  the  real  axis  are 
easily  obtained.  Define  values  of  the  complex  variable  s  on  the  real  axis 
by 

sk  -<T  (14) 

For  values  of  s  on  the  real  axis  the  characteristic  equation  (4)  becomes 

f(o<,/S,  <T)  =2  (**>*+■  £Ck  +  e0 (--«")*  (15) 

K=o 

Of  particular  importance  is  the  degeneration  of  the  characteristic  equation 
on  the  real  axis  into  a  single  real  equation  whose  imaginary  part  is  identi- 
cally zero.  For  easier  discussion  define: 

N  ic 


B_((T)  -T("l)  L   (r1 

CR((r)  =  X/"1^0*  ^  (15a) 


fc=o 

k=o 
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and  write  equation   (15)    as 

o<6R(r)-h/?CR(r)  -hE*R(<r)  =  o  <ie) 

Equation  (16)  is  a  single  equation  in  three  unknowns  o(  ,  j3  and  5"  which  on 
selection  of  a  real  point  ((T"0  )  to  be  mapped  defines  a  curve  in  the  o(-^ 
plane.  For  the  case  under  discussion  o(  and  &   appear  linearly  and  the 
curve  will  be  a  straight  line. 

Equations  (13)  and  (16)  are  the  keys  which  allow  interpretation  of  the 
characteristic  equation  in  the  parameter  plane.   Equation  (13)  is  applicable 
for  all  s-plane  points  except  those  on  the  real  axis  and  it  maps  points  into 
points  and  s-plane  contours  into  o(-j8  plane  contours.   Equation  (16)  represents 
the  degenerate  case  when  s  assumes  a  value  on  the  real  axis  and  it  maps  a 
single  s-plane  point  into  a  straight  line  in  the  o(-|?  plane.  The  significance 
and  interpretation  of  these  parameter  curves  is  explained  in  detail  in  the 
next  chapter. 

Parameter  Plane  Equations  (  o<^  Product  Case) 

All  too  often  system  parameters  appear  in  the  characteristic  equation 
non-linearly;  an  occurence  which  precludes  parameter  plane  analysis  by  the 
techniques  thus  far  developed.  Theoretically  parameter  plane  techniques  can 
be  extended  to  cover  any  eventuality  but  the  impracticability  if  not  the  im- 
possibility of  obtaining  explicit  solutions  for  system  parameters  in  terms 
of  s-plane  root  locations  prevents  theory  from  becoming  practice.  Character- 
istic equations  in  which  ©</?  products  appear  are  a  frequently  occuring  and 
important  exception  to  this  statement. 

f(o<  tjS  ,  S  )  -^  (*bK  +  /6cK  +  o^ci<  +  eK)SK  =  O  (17) 


K=0 


13 


Parameter  plane  notation  and  definitions  for  the  oUS     product  will  be  consist- 
ent with  those  employed  for  the  linear  case.   Noteworthy  is  the  continuing 
relevance  of  relations  (5)  through  (9).   Separation  of  equation  (17)  into 
its  real  and  imaginary  parts  yields: 

f  (-<)'  («>>*+  /8cK  +  -pVek)  H"  T*  (ft)  -  o 

MO 

J  (-*)"  (°<  tR  +  (8  cK  +  x^d,  +  e» )  h"  Uk  (ft)  =  o 

IfcO 

The  same  manipulations  employed  in  the  linear  case  allow  these  equations  to 
be  expressed  as: 

J  (-»)"  («*>k+  /8cs+  o^JK  +  e«)  ^;  UM  (S)  -o 


Y 


In  the  interests  of  brevity  define: 

»1<s.^-2(-0\k''Uk-,(s) 

K*o 


b2(  £  .  H, 
cx(  5  .  H* 
c2(  5  -  uM 

V  $  .  H, 
e2<£  ,w„ 


(-»)  c,utfU.(ft) 

JtsO 

-£(-')*  d^lW  ft) 

-t  (-0*4^^  (5) 

K=0 

-  f  (-0*  ek  M?  U,  (ft) 


n=o 


and  rewrite  equations   (18)   as: 
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(18) 


(19) 


(20) 


(Notice  that  definitions  (19)  are  identical  to  definitions  (12)  with  D  (  £  , 
LJN}   and  D  (£  ,  Uu   )  added).   Equations  (20)  can  be  solved  for  o(  and  /S 
as  functions  of  £  and  (jj      by  three  different  methods. 

1.  The  parameter  /S  is  initially  eliminated  from  equations  (20) 
and  0<  is  obtained  from  the  single  resulting  equation  using 
the  quadratic  formula.  Values  for  /Q  are  subsequently  found 
from  a  linear  formula  utilizing  the  previously  found  values 
of  0<  as  knowns.   Starting  with  equations  (20) 

O^CB^  -  B2D1)  +<*(B1C2  +  D2E1  -  B2C1  -  D^)  +  (C^  -   C^)  =  0 

'  ^ftCb-flkO,))  (21) 

-   teeC.  +  qEp-B.Ca-OaE;)5-  V  (8,  Dt-£g  flJfcaF.  -C,F») 

^CftH)-— ^Jf-        ^^H)=-^^      (22) 
The  last  two  formulas  yield  identical  results  for  ft  . 

2.  This  method  duplicates  method  1,  except  that  0<  is  initially 
eliminated  and  yQ  is  solved  via  the  quadratic  formula. 

#(;M)  =  (8,C^D,E,-BgC,-D.E,)     +  («) 

^(QDe-CgO,; 
±l/(8£C  +  DeF,-6)Ca-Dl£'g)g  -  4  (C.Pg-CaDjfSefi-aggT 

~  (24) 

3.  A  third  procedure  involves  solution  for  both  c<(  £  »  UJN   )  and 

j8  {  C  9  U)N   )  by  the  quadratic  formula  after  3    and  <?(.  , 
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respectively  are  eliminated  among  equations  (20).  The  re- 
sults are  solutions  (21)  and  (23)  which  will  be  designated 

c*(£  ,  UH)  =  X±{T 
/8<£  9UJH)  «  U±fT 

While  equations  (21)  through  (24)  display  explicit  solutions  for  oC 
and  /S    a  relevant  question  still  remains  -  are  solutions  (21)  and  (22)  iden- 
tical to  those  obtained  in  (23)  and  (24)?  This  question  paraphrased  to 
apply  to  the  third  procedure  would  be  -  are  each  of  the  possible  point  sets 
( ©(  9   S   )  valid  solutions  for  equations  (20)  or  are  two  of  these  sets  ex- 
traneous? Direct  substitution  reveals  that  the  solutions  obtained  using 
equations  (21)  and  (22)  are  identical  to  the  solutions  obtained  using  equa- 
tions (23)  and  (24).  Moreover,  using  the  third  procedure  yields  valid  solu- 
tions: 

( 0(  ,  |S  >  -  (x  +  (T ,     u  -  /▼)   and  (  0<  ,  (3  )  =  (x  -  fT,  u  +  V^) 

The  two  remaining  possibilities  are  extraneous  and  will  not  satisfy  equations 
(20). 

As  for  the  linear  case  s-plane  points  on  the  real  axis  are  degenerate; 
the  imaginary  restriction  on  the  characteristic  equation  is  identically 
satisfied  and  only  one  equation  relates  the  three  variables  o(t   /3    and  s. 
This  equation  which  is  appropriate  for  points  on  the  real  axis  is  obtained 
by  defining  s  as  in  (14)  and  substituting  into  characteristic  equation 
(17). 

f(<*  y  (3  ,o~)  «]r  (*/>K^ck+oy3dK-he<)  (-<r)   =o    <25> 

K=0 

As  before  define: 

R 


c5(r)  -y  (-\)kcK<r*  (25a) 


/Y  U 

dr((d-  y  (-0  dK°"K 

and  write  equation  (25) 


as: 


/3  (*,  (T)  =  _  «B*«T)  -hER(<T)  (26) 

The  locus  of  this  equation  for  any  given  value  of  ft~*  describes  a  hyperbola 
in  the  parameter  plane  with  asymptotes  of: 

o<  =  _  C«(r)  and         a  „  _  B„(<r) 

DR«r)  DR(c-) 

Comparing  these  explicit  expressions  with  those  occuring  for  the  case 
of  coefficients  linear  inland  fS   we  find  three  mentionable  differences. 

1.  In  the  linear  case  all  s-plane  values  mapped  into  real  values 
in  the  o(-f3  plane.   For  the  o0#  product  case  this  is  no 
longer  true;  the  possibility  for  complex  0(,/?  values  exists. 
Physically  this  implies  that  it  is  impossible  to  place  s- 
plane  roots  at  such  locations  by  adjustment  of  pC  and  A3. 

2.  In  the  linear  case  each  complex  s-plane  point  mapped  into  a 
single  point  in  the  parameter  plane.  For  oOS  products  a 
complex  s-plane  point  usually  maps  into  two  different  points 
on  the  vL-3  plane. 

3.  For  both  linear  and  product  cases  real  s-plane  points  have 
curves  as  <*>-0  plane  images;  however,  where  the  curve  is  a 
straight  line  for  the  linear  case  it  becomes  a  hyperbola  for 
the  product  case. 

Parameter  Plane  Methods  Extended  Beyond  the  Product  Case 

Extension  of  parameter  plane  methods  beyond  the  product  case  appears 
prohibitive  even  when  explicit  solutions  for  o(  and  f3     can  be  found. 
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Suppose  we  allow  characteristic  equations  of  the  form 

k=o 
Under  these  assumptions  explicit  solutions  for  o(  and  &     exist  and  involve 
no  insurmountable  or  extraordinary  difficulties;  in  fact  simultaneous  solu- 
tion of  two  conies  for  o(   and  @     yields  the  desired  expressions. 

**A11  +  «P  A21  +    (*\l   +    *A41  +    Z8  A51  +   A61  =  ° 

°^A12  +    *£  A22  +  ^A32  +         *  A42  +         P  A52  +       A62  =  ° 

where 


An<£.wN)-^  (-')  CtIKUM  UK_,  (5) 
a32(  c, ,  w„  )  -  ^  (-/)"  Q3K  uiZ  UK  (5) 


KsO 

The  problem  of  consequence  is  interpretation  of  the  parameter  plane  curves 
since  each  s-plane  point  generally  has  four.o<-/3  plane  images.   Additionally, 
the  four  valid  solutions  for  ( oC  >  £»  )  must  be  distinguished  from  among  a 
minimum  of  eight  generated  in  the  process  of  solving  the  simultaneous  equa- 
tions. 

Parameter  Plane  Methods  for  More  Than  Two  Parameters 

The  parameter  plane  is  a  competent  two  variable  technique  in  analogy  to 
the  root  locus  being  a  capable  single  variable  technique.   Pursuing  this 
analogy  we  find  that  parameter  plane  technique  can  accommodate  three  variables 
with  a  facility  about  equal  to  that  of  the  root  locus  to  accommodate  two 
variables.   Problems  having  a  third  parameter  plane  variable  can  be  approach- 
ed in  three  ways. 

1.   An  additional  restriction  or  specification  is  invoked  (such 
as  a  steady  state  error  restriction)  which  permits  solution 
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for  the  third  variable,  ^f  ,  as  a  function  of  the  remaining 
two,  f(  oC  ,  jo  )•  With  elimination  of  the  third  parameter  the 
problem  reverts  to  the  standard  two  parameter  problem  and 
analysis  proceeds  as  before.  However  linearity  of  oC  and  /S 
in  the  coefficients  of  the  characteristic  equation  is 
commonly  destroyed  when  eliminating   }f   in  favor  of  f (  0<  , 
/3  ).   This  shortcoming  seriously  diminishes  this  otherwise 
attractive  approach  since  loss  of  linearity  usually  precludes 
parameter  plane  analysis. 

2.  Addition  of  a  third  dimension  to  the  two  dimensional  para- 
meter plane  allows  adequate  representation  of  the  three 
variable  problem.  Under  this  representation  an  s-plane 
contour  maps  into  a  surface  in  (  oC  »  /Q  >  0  )  parameter 
space.  Any  choice  of  system  parameters  which  lie  on  this 
surface  will  have  at  least  one  pair  of  roots  on  the  chosen 
s-plane  contour.  Only  the  problems  innate  to  three  dimensional 
presentation  limit  the  attractiveness  of  this  method. 

3.  One  selected  s-plane  contour  can  be  mapped  into  several  o(.-/3 
plane  curves,  each  curve  representing  a  different  value  of  tf  , 
the  third  parameter.  This  procedure  effectively  extends  para- 
meter plane  analysis  to  three  variables  especially  when  only 
one  or  two  s-plane  contours  are  of  interest.  Another  credit 
is  that  applicability  is  not  limited  by  the  manner  in  which 
the  third  parameter  appears  in  the  characteristic  equation. 
Moreover,  the  o(-^w  plane  image  of  any  specified  s-plane  point 
(real  axis  points  excluded)  plots  on  a  straight  line  as  the 
third  parameter  is  varied  whenever  all  parameters  are  linear 
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with  respect  to  one  another. 
This  last  observation  is  quite  useful  in  the  following  context.   Suppose  we 
specify  any  s-plane  point  as  s(x,y)  where  x  is  held  constant  and  y  is  the 
running  variable  along  the  s-plane  contour  of  interest.   The  map  of  this 
contour  in  the  oC-^S  plane  depends  on  the  third  parameter;  a  different  mapping 
for  each  different  value  of  o  .  Connect  the  points  which,  on  each  curve, 
represent  the  same  value  of  the  running  parameter  y  -  the  locus  of  these 
points  is  a  straight  line  in  the  parameter  plane.   Furthermore,  the  spacing 
of  the  various  image  curves  is  proportional  to  the  value  of  0  . 

The  apparent  utility  of  these  linearity  properties  warrants  the  more 
rigorous  foundation  supplied  by  formal  theorem  and  proof. 

Theorem  I.   Suppose  three  variable  parameters  appear  linearly  in  the 

coefficients  of  the  characteristic  equation. 

f(oC  ,£,*,■)  =2^(°<tK^/3Ck-hy  ^k  +  (2k)SK  =  0     (27) 

K=o 
Then  the  locus  of  any  fixed  s-plane  point  in  the  parameter  (  &£ — B   )  plane 

is  a  straight  line  and  the  image  of  the  point  is  a  linear  function  of  the 

third  parameter 

Proof.   Resolve  equation  (27)  into  its  real  and  imaginary  parts,  equate 

each  part  to  zero  and  simplify  as  was  done  when  originally  formulating  the 

parameter  plane  equations.   The  result  is  two  equations  analogous  to  equations 

(11). 

oCB^  £,  WN)  +  pcx(  £  ,(JN)  +  Yf1(  £  ,  UH)  +  E1(  J  ,UJH  )  =  0 

(28) 
<*B2(  5   ,  U/N  )  +  /3B2(  £   ,  U)M  )  +  YF2(   £  ,  Oi  )  +  E2(  £   ,  LUN  )   =  0 

where 

V  $  •  wN )  -  J  (~ i) "  *k  U>*  U,_,  ( £)  j 

K  =  0 
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(2(^)=J(-/)KfKu;wKUK(c) 


and  all  other  quantities  are  as  defined  in  (12) 

Explicit  solution  of  equations  (28)  yields  the  customary  parameter  plane  equa- 
tions. 

(29) 


For  a  fixed  point  in  the  s -plane  £  and  ldN  become  knowns  and  cL  and  ^be- 
come functions  of  the  single  variable  Y  .  Under  these  conditions  equations 
(29)  are  recognizable  as  the  parametric  representation  of  a  straight  line  in 
the  oC-f2  plane.  Furthermore  the  linear  appearance  of  )f  in  both  of  these 
explicit  expressions  implies  that  the  position  of  the  image  point  on  the 
straight  line  locus  is  linearly  proportional  to  the  value  of  Y  .   Q.E.D. 

Parameter  Plane  Equations  in  the  Z-Plane 

The  aforementioned  analogy  between  parameter  plane  methods  in  the  s  and 
z-planes  can  be  strengthened  to  practical  identity  of  procedures  and  equa- 
tions.  (The  only  noteworthy  difference  between  the  two  planes  is  choice  of 
contours  to  be  mapped  onto  the  o<-^  plane;  constant  £  contours  are  popular 
in  the  s-plane  whereas  spirals  of  constant  damping  or  circles  of  constant 
settling  time  are  popular  in  the  z-plane).   This  close  correspondence  made 
formulation  of  the  parameter  plane  equations  in  the  more  familiar  s-plane 
with  subsequent  extension  to  the  z-plane  the  most  attractive  alternative. 
For  z-plane  formulation  we  assume  a  characteristic  equation  of  the  form 

f(oc./3. «)  =^(°itl|(+/scK-KeK)zK  =o 

K=o 
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with  provision  for  inclusion  of  terms  for  the  oC/6  product  case,  three 
parameter  case,  etc.   Results  derived  in  the  s-plane  are  valid  in  the  z-plane 
subject  to  the  following  modifications. 

a)   If  the  complex  variable  z  is  defined  as 


b) 


•UJj-^+df^) 


With  z  so  defined  no  modification  of  s-plane  formulas  is 
necessary  beyond  replacing  4  by  £   and  U)N   by  lu^  . 
If  the  complex  variable  z  is  defined  as 

•■'^(ft.  +  Jy/'-ft*    ) 

The  total  effect  of  so  defining  z  is  to  obviate  the  need  for 

1/ 

the  (-1)    factor  in  the  definitions  of  B,C,D,E  and  F. 

Specifically,  all  s-plane  formulas  are  valid  if  the  definitions 
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9)    (25a)   and   (28)   are  redefined  as 
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The  reader  is  cautioned  that  the  analogy  between  C   and  £  and  U^  and  U^ 
extends  only  to  the  mathematical  formulation  of  the  parameter  plane  equations. 
The  value  of  £_  is  in  no  way  indicative  of  system  overshoot  nor  is  Cl/. 
indicative  of  the  system's  undamped  natural  frequency. 

Mitrovic's  Method  and  Extensions  of  Mitrovic's  Method 

Although  Mitrovic  methods  are  special  cases  of  the  parameter  plane  they 
warrant  mention  because  of  a  relation  between  the  o(  -  (3     (M-plane)  curves  and 
the  s -plane  roots  not  present  in  the  more  general  case.  This  relation  and 
its  application  to  the  z-plane  will  be  discussed  in  Chapter  3;  the  present 
objective  is  simply  to  define  Mitrovic  methods.   Parameter  plane  theory  assumes 
a  characteristic  equation  of  the  form: 

N 

f<*,/8,  s  )  -£(otbk+/9cK  -heK)sK  =  o 

Mitrovic's  method  and  its  extensions  assume  that  only  one  b.  and  one  c. 
(i  4   j)  have  non-zero  value;  oC  and  &    are  subsequently  redefined  as: 
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(3'-  /S  Cj 

The  characteristic  equation  in  a  format  apropos  to  Mitrovic  formulation  then 


appears  as: 

K=0 
If  i  ■  0,  j  =  1  the  formulation  is  the  one  originally  advocated  by  Mitrovic; 


;.,,'..>.£  eKsK  +    *V'mV-0 


otherwise  the  formulation  is  referred  to  as  an  extension  of  Mitrovic 's  method, 
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CHAPTER  III 

The  Interpretation  Problem 

To  regard  a  functional  relation  between  two  complex  variables,  w  ■  g(z), 
as  limited  to  defining  w  in  terms  of  z  is  indeed  tunnel  vision.   Lost  would 
be  such  concepts  as  images  of  z-plane  contours  and  mappings  of  z-plane  areas 
into  selected  areas  of  the  w-plane.  An  equally  grievous  error  would  be  to 
regard  the  explicit  expressions  for  oC(   £  ,  UJM  )  and  XJ  (  £  ,  UJN  )  as  formulas 
to  place  roots  at  specified  s-plane  (z-plane)  locations.  Such  regard  would 
obscure  two  most  important  interpretations,  that  of  mapping  and  that  of  parti- 
tion of  the  oC_ 0    plane  into  regions  having  a  similar  root  property.  The 
problem,  therefore,  is  not  for  a  single  interpretation  of  oL-B  plane  curves 
but  rather  it  is  for  an  answer  to  the  question,  "How  can  the  significant  in- 
formation buried  in  s-plane  (z-plane)  pole  zero  configurations  be  extracted 
from  parameter  plane  curves?" 

Interpretation  in  the  &--fi   Plane 

A  simple  yet  comprehensive  answer  to  the  interpretation  problem  is  pro- 
vided by  examining  various  point  sets  in  the  o(r0   plane  for  common  properties. 
In  figure  (la)  a  z-plane  contour  terminating  at  two  points  on  the  real  z- 
axis  is  designated;  figure  (lb)  shows  the  oC-8  plane  image  of  this  contour 
and  the  images  of  the  two  real  z-plane  points  terminating  this  contour.  With 
reference  to  figure  (lb)  define  the  following  point  sets. 

Point  Set  A  -  Point  set  A  contains  the  images  of  all  points  on  the 
designated  z-plane  contour. 

Note:  This  set  contains  all  points  on  the  image  curve 
drawn  in  figure  (lb)  and  no  points  not  on  this  curve. 
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Point  Set  B  -  Point  set  B  contains  all  points  which  are  the  image  of 

the  real  z-plane  points  5Z    and  6Z   . 
Point  Set  C  -  Point  set  C  contains  all  points  in  the  cross-hatched  area 

and  it  is  bounded  by  but  disjoint  from  set  A. 
Point  Set  D  -  Point  set  D  contains  all  points  in  the  dashed  shaded  area 

and  it  is  disjoint  from  set  A  and  set  B. 


z-plane 


closing 

conjugate 

contour 


contour 


FIGURE  la 


FIGURE  lb 


Let  M(  oC  ,  /3   )  be  any  point  in  the  parameter  plane  and  consider  the  signi- 
ficance of  M  belonging  to  certain  point  sets. 

Observation:   If  M  belongs  to  set  A  the  characteristic  equation  must 
have  one  complex  root  lying  on  the  chosen  z-plane  contour  and  another  root 
lying  on  its  reflection  in  the  real  axis.   Recall  that  explicit  equations 
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give  the  parameter  values  needed  to  place  a  root  at  (  £   ,  U)N     ) . 

However,  these  explicit  equations  can  also  be  considered  as  mapping  functions 

which  establish  a  relation  between  z-plane  points  and  oC-@   points.   Regarded 

thus,  every  point  in  set  A  has  its  domain  on  the  designated  z-plane  contour 

and  satisfies  these  explicit  expressions;  hence  the  conclusion  observed  is 

inevitable. 

Observation:   If  M  is  contained  in  set  B  the  characteristic  equation  must 
have  at  least  one  real  root  at  z  =  (JZ      (or  (Z       )  . 

Observation:  All  M-points  located  in  set  C  have  an  equal  number  of 
characteristic  equation  roots  enclosed  by  the  z-plane  contour  and  its  reflec- 
tion in  the  real  axis.   Furthermore  the  number  of  characteristic  equation 
roots  enclosed  by  an  M-point  in  set  C  differs  by  two  from  the  number  of  roots 
enclosed  by  an  M-point  in  set  D.  The  validity  of  this  statement  can  be 
established  by  contradiction.  Assume  two  points  exist  in  set  C  which  have  a 
different  number  of  characteristic  equation  roots  enclosed  by  the  z-plane 
contour.   Then,  by  continuous  motion  of  the  M-point  from  one  of  these  points 
to  the  other,  characteristic  equation  roots  could  move  across  the  closed  z- 
plane  contour  without  ever  once  having  the  roots  lie  on  the  contour  itself. 
(All  M-points  having  roots  on  the  contour  belong  to  set  A  which  is  disjoint 
from  set  C.)  This  is  impossible;  therefore  by  contradiction  all  points  in  C 
have  the  same  number  of  characteristic  equation  roots  enclosed  by  the  z-plane 
contour.   Similar  continuity  arguments  can  be  applied  to  show  that  a  differ- 
ent number  of  roots  are  enclosed  for  M-points  in  set  C  and  set  D.   (To  deny 
this  proposition  would  imply  that  a  root  lying  on  the  specified  z-contour  could 
be  moved  in  any  direction  in  the  z-plane  without  affecting  the  number  of  roots 
enclosed  by  the  contour).  That  the  enclosed  roots  in  set  C  and  D  differ  by 
two  stems  from  the  restriction  that  coefficients  of  the  characteristic  equa- 
tion be  real;  root  motion  across  the  contour's  upper  limb  requires  motion  of 
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the  conjugate  root  across  the  contour's  lower  limb. 

Observation:  The  &/.-&    plane  images  of  points  at  which  the  z-plane 
contour  intersects  the  real  z-axis  also  serve  to  divide  the  parameter  plane 
into  point  sets.  M-points  belonging  to  sets  on  opposite  sides  of  such  con- 
tours will  differ  by  one  in  the  number  of  characteristic  equation  roots  en- 
closed by  the  z-plane  contour.   Justification  of  this  statement  is  analogous 
to  that  of  the  preceding  observation. 

Preceding  observations  although  they  cannot  claim  mathematical  rigor 
provide  a  framework  for  parameter  plane  interpretation  one  step  above  the 
intuitive  approach  usually  advanced.  For  application  these  observations 
can  be  profitably  restated. 

1.  M-points  on  an  image  curve  imply  a  characteristic  equation 
root  on  the  z-plane  (s-plane)  contour  upon  whose  image  the 
M-point  lies 

2.  Parameter  plane  images  of  z-plane  contours  and  images  of 
the  intersection  points  of  these  contours  with  the  real 
axis  serve  to  divide  the  <^— 8   plane  into  regions  or  point 
sets.  The  distinguishing  characteristic  of  each  region  is 
that  all  M-points  within  a  region  have  the  same  number  of 
characteristic  equation  roots  enclosed  by  the  z-plane  (s- 
plane)  contour. 

3.  If  the  boundary  between  regions  is  the  image  of  a  complex 
contour  the  number  of  roots  enclosed  for  M-points  located 

in  adjacent  regions  differs  by  two.   If  the  boundary  between 
regions  is  the  image  of  a  real  axis  point  the  number  of  roots 
enclosed  for  M-points  in  adjacent  regions  differs  by  one. 
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Shading  Convention  and  Rules 

The  interpretation  of  d--f?   curves  as  images  of  selected  s-plane  (z- 
plane)  countours  was  previously  discussed.  Moreover,  the  preceding  para- 
graph showed  that  movement  of  the  M-point  across  an  oL-fi   curve  corres- 
ponds to  two  complex  roots  moving  across  the  selected  s-plane  (z-plane) 
contour.   The  direction  in  which  the  roots  move  across  the  s  plane  contour 
is  usually  available  by  inspection  of  the  direction  in  which  the  M-point 
moves  across  the  oC-y8  curve.   (i.e.  ,  Motion  of  the  M-point  across  the  o(-/8 
plane  image  of  the  contour  £  =  £  in  the  direction  of  oi.-/9  plane  images  of 
contours  £  >.  £tf  implies  that  the  s-plane  roots  move  across  the  contour  in 
the  direction  of  increasing  £  ) .  However,  parameter  plane  regions  in  which 
images  of  two  different  C  contours  intersect  or  in  which  the  image  of  a 
given  s-plane  contour  loops  on  itself  do  not  readily  allow  the  direction 
of  root  motion  across  the  s-plane  contour  to  be  inferred  from  the  motion  of 
the  M-point  across  its  oC-fi  plane  image.   For  such  cases,  and  in  fact  when- 
ever required,  the  direction  of  root  motion  across  an  s-plane  contour  as 
the  M-point  moves  across  its  image  in  the  parameter  plane  can  be  deter- 
mined mathematically. 

Adopt  the  convention  that  the  parameter  plane  images  of  constant  £ 
(or  constant  kJN     )  s-plane  contours  are  shaded  so  that  an  M-point  cross- 
ing the  oC-0    curve  from  shaded  to  unshaded  side  causes  two  complex  s-plane 
roots  to  leave  that  portion  of  the  s-plane  enclosed  by  the  contour. 

Theorem  I  The  sign  of  the  Jacobian   T{   ' '   J  determines  the  side 
of  the  image  curve  to  be  shaded.   If  the  Jacobian  is  positive  (negative) 
the  left  (right)  hand  side  of  the  image  curve  should  be  shaded  facing  in 
the  direction  of  increasing  UJu  (  £  ). 

Proof   Construct  a  cartesian  three  dimensional  vector  space  on  the 
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parameter  plane  whose  position  vector  is: 

r  (  oi. ,  {8 ,  2  )  -  oCx-h/Oj-r-zlT 

In  this  space  all  curves  lying  in  the  parameter  plane,  as  all  images  of  s- 
plane  contours  must,  will  have  position  vectors  of  the  form 

r<ot,/G)  =  o(.X+  fit 

Now  recall  that  a  requirement  for  the  construction  of  parameter  plane  curves 
was  the  existence  of  explicit  solutions  for  oC  and  |8  in  terms  of   C  and 
U)N      .   Namely: 

This  requirement  allows  the  position  vector  r(  eC  ,  /S    )  to  be  written  as 

r<£  ,14 )  -  oc<  $  ,  H,)/  +   /*  (  £  .<*t>T 

The  oC-/3  plane  image  of  a  constant  j?  contour  of  value  :?«>  is  traced  by 
the  vector 

Fix  attention  at  one  point  on  this  curve  and  determine  the  motion  of  that 
point  resulting  from  an  incremental  increase  in  £>   .It  is 

Obviously  the  direction  of  M-point  motion  in  the  direction  of  increasing  > 
at  this  point  is  determined  by  the  vector  2-Llatj  zSL.      Since  an  M-point 
moving  in  a  direction  to  increase  £  causes  roots  to  enter  the  area  en- 
closed by  a  constant  £   contour  that  side  of  the  oC—B     plane  image  toward 
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which  the  vector 


*r~  points  1 


s  shaded.   Graphically  this  scheme 


suffices  but  it  can  also  be  stated  in  a  mathematical  formulation.   By 
definition  the  tangent  to  the  o<,-  S    plane  image  of  the  constant   £=>o  s~ 
plane  contour  at  the  point  (  >ft  ,  i*J0    )  in  the  direction  of  increasing 
U        is  given  by  the  vector  SlCLggii^gJ  .   If  the  counter-clockwise  rotat 
of  the  vector  °^r^i^  into  the  vector  °r(  \°\     *'    is  more  than  180°  the 
latter  vector  points  to  the  right  of  the  constant  £0  curve  image  and  that 


ion 


side  of  the  o£-^6  curve  should  be  shaded.   Contrariwise,  if  the  rotation  is 
less  than  180°  the  vector  ^r'^°  >  °'   points  to  the  left  of  the  constant  £0 
mapping  in  the  oL—p     plane  and  that  side  of  the  curve  should  be  shaded. 
Fortunately  this  criterion  is  the  exact  definition  of  the  cross-product  of 
two  vectors  lying  in  a  cartesian  plane. 

dm        A  71 

If  A  is  negative  shade  the  right  hand  side  of  the  oL— 8    plane  curve  facing 
in  the  direction  of  increasing  UJN   ;  if  it  is  positive  shade  the  left  hand 
side.   At  this  point  reconcile  the  sign  of  A  with  a  more  familiar  expression. 
To  do  this  express  r(  C,  ,  UJ    )  in  terms  of  its  <>C  and  &   components  and 
take  the  required  derivatives. 

Jr(C0,iuO  .  Joc(g..t4,)  ;f  ,      J/3 te, h>0  "J 


ete 


and 


Jr(£oluJo)  y   Jr($otWo)     m 
3uJH       *  ()$ 


c)uJ„ 


£$ 
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Hence  the  sign  of  A  depends  on  the  sign  of  the  determinant  which  is  by 
definition  the  well  known  Jacobian  ~T(  =-5). 

Extension  of  these  results  to  other  s-plane  contours  and  also  to  z- 
plane  contours  is  possible.   In  fact  some  general  theorems  must  exist 
which  encompasses  all  alternatives,  but  the  statement  of  such  a  theorem  and 
the  determination  of  necessary  and  sufficient  conditions  are  difficult. 
Since  there  are  but  a  limited  number  of  commonly  used  s-plane  (z-plane) 
contours  it  suffices  to  list  the  governing  Jacobian  for  each.   The  proof  in 
each  case  proceeds  exactly  as  that  presented  above. 

Extension  of  Theorem  I    The  sign  of  the  Jacobian   7"  (  y  y  )  >  where 
the  independent  variables  x,y  are  dependent  on  the  contour  being  mapped, 
determines  the  side  of  the  image  curve  in  the  ot-(3    plane  to  be  shaded. 
Shading  rules  are  as  listed  in  theorem  I. 

a)   For  constant  settling  time  or  frequency  contours  in  the  s- 


Ati) 


plane  the  governing  Jacobian  is 


'H 


dot    JuJ„  ¥  c)*  £_£_ 


where  s  =  U)„  (-  £  -h>j^j-£2  )=  -  $-   +■  $uJ 
b)   For  constant  £  or  \jj      (settling  time)  contours  in  the  z- 


plane  the  governing  Jacobian  is 


Jfl 


where  *  -  UJ     (- C,  +  ^1-$*    ) 
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3ft$)- 


Note:   Frequently  z  is  expressed  as  z  *»  (//.  (Q    -\-  $  y]  -  C,     ) 

in  which  case  the  sign  of  the  Jacobisn  is  reversed. 

c)   For  constant  damping  contours  in  the  z-plane  the  governing 

Jacob ian  is 

dot      gljdto  JoC_  ijb  ofc      o)^,       J/3     c)$2 


doc    Juki   .  ^oC  o>C-, 


Mathematical  formulation  for  shading  is  virtually  complete  except  for 
the  degenerate  case  of  points  on  the  real  axis.   Each  of  these  points  have 
an  entire  curve  as  their  image  in  the  parameter  plane  since  only  one  equa- 
tion exists  (16)  which  relates  the  three  variables  oC,  8  ,  S"  •   (The  restric- 
tion that  the  imaginary  part  of  the  characteristic  equation  equal  zero  is 
identically  satisfied  by  points  on  the  real  axis) .   If  parameter  plane  curves 
can  be  constructed  explicit  solution  of  equation  (16)  for  oC  in  terms  of  /& 
and  (Tor  for  j8    in  terms  of  0(.  and  (T  exists.   Designate  these  explicit 
expressions  as 

<*=  <*(£,$") 
/3  mfi  (oc)(r) 
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Theorem  II   The  shading  of  the  parameter  plane  curve  which  is  the 
image  of  a  real  point  (J~(  0"k  )  in  the  s-plane  (z-plane)  is  determined  by 
the  sign  of  the  partial  derivative  ^fi  (*o  >  O  (or  ^gLLfjgj  G~J   ^   where 

(  o£0  ,  R      )  is  a  point  on  the  parameter  plane  image  of  (yQ        .   The  shading 
rule  is  as  follows. 

(a)   If  increasing  (T  (  or  (J^   )  from   0^   causes  a  real  root 
to  enter  the  area  enclosed  by  the  contour  shade  the  left 


(right)  side  of  the  curve  facing  in  the  direction  of  increas- 
mg   OC   when  -^ — ^ — '  ■■■- 
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is  positive  (negative). 


Equivalent  to  this  is  the  statement  that  if  increasing  C* 
causes  a  real  root  to  enter  the  area  enclosed  by  the  contour, 
shade  the  left  (right)  hand  side  of  the  curve  facing  in  the 

a-  .•      c     •        •     />     u  oV.(&.<To) 

direction  of  increasing  S   when  - =-* — « — -  is  negative 

1  o?<T 

(positive) . 
(b)   If  increasing  (T  (or  6"^  )  from   6"^   causes  a  real  root  to 
leave  the  area  enclosed  by  the  s-plane  (or  z-plane)  contour 
shade  the  right  (left)  side  of  the  oL-(S    curve  facing  in  the 

j-  r  •  «/    u     Jfi   (<*■<?  >  (To)   • 

direction  of  increasing  Oc  when  — ' — ~ — ' is  positive 

(negative) . 
For  comparison  a  table  presents  this  theorem  in  a  most  effective  manner. 


Case 

face 
Direcrion 

S/gn    of 

0)0- 

Shade 

face 
Direcrion 

Siqn    of 

Shade 

o)(T 

Increased    <T 
(or  aj)    implies 
Cl   coot   enrers 
f/>t    enclosed 
area 

Increas- 
ing 

Po  s  1  ft  v  e 

Left 

Increas- 
ing 

A/egaftve 

Left 

A/egaftve 

Rignt 

Positive. 

/?ignt 

Increased    (T 
(or  (f)    implies 
a.  root  leaves 

iht.    enclosed 
area 

Increas- 
ing 

Positive 

Right 

Increas- 
ing 

fi 

A/egat/i/e 

Right 

A/egaftve 

Left 

Pos/ftve 

Left- 

Table  I 
34 


Values  on  the  real  s-axis  and  the  real  z-axis  are  not  defined  identically; 
the  usual  case  being 

Real  s  ■  —  <T* 

Real  z  =  4-  (T^ 
The  manner  of  definition  does  not  affect  the  validity  of  the  theorem  if  the 
user  is  consistent  throughout.   For  instance,  under  these  definitions  a 
change  of  s  from  -1  to  -2  would  indicate  an  increase  in  (T*  whereas  a 
similar  change  in  z  would  indicate  as  a  decrease  in  (J"-g.    .   Notice  also 
that  either  — — - — ' — 6  or  Hr— ' fc  carry  sufficient  information  to 

cP<r  o><r 

shade  and  their  dictates  must  always  agree.  An  exception  occurs  when  the 
tangent  to  the  oL-/S    curve  is  vertical  or  horizontal  in  which  case  one  of 
the  partials  will  degenerate  to  a  value  of  zero  or  infinity  while  the  other 
will  still  impart  useful  information. 

The  proof  to  this  theorem  is  analogous  to  that  presented  for  theorem  I. 
In  this  case  either  oC      and   fiT*  or    /3         and  (f     can  be  considered 
independent  variables  depending  on  which  variable  (  OL       or   |S   )  is 
explicitly  expressed.   The  position  vectors  are 

P(oC,(r,)=sO^/.  +  /G(0N*>0OJ  for  explicit  expression  of  ft 

— k  ^  -*-  — *■  i 

P  (  p  »  <T*  )  =   oc(|S>0")/<.+  /Q^   for  explicit  expression  of  06 

From  this  point  forward  the  proof  is  identical  with  that  of  theorem  I. 

A  few  concluding  comments  on  shading  and  shading  formulas  should  be 
reiterated  for  the  reader.   First,  the  shading  of  parameter  plane  curves 
provides  a  comprehensive  picture  of  root  motion  across  selected  s-plane 
contours  as  a  function  of  the  simultaneous  variation  of  two  system  para- 
meters. When  several  s-plane  contours  are  mapped  onto  the  o£-/2  plane  this 
root  motion  is  usually  evident  without  the  tedium  of  shading  and  if  shading 
is  desired  it  can  be  accomplished  by  inspection.   For  any  case  where  inspec- 
tion shading  is  impossible  the  formulas  presented  here  resolve  the  difficulty. 
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In  this  they  extend  beyond  existing  criteria  which  can  be  applied  only  for 
coefficients  sl    linear  in  oC  and  /Q       [4],   However,  the  labor  involved 
in  evaluating  the  sign  of  a  Jacobian  is  somewhat  greater  than  that  required 
for  existing  shading  procedures  and  the  former  undoubtedly  will  find  employ- 
ment only  as  a  last  resort.   Probably  the  most  significant  contribution  of 
these  rules  is  that  they  put  shading  on  a  more  rigorous  framework. 

Example  of  Parameter  Plane  Curve  Interpretation. 

An  example  will  serve  to  illustrate  the  salient  features  of  parameter 
plane  interpretation  and  curve  shading  discussed  in  the  preceding  sections. 
Consider  the  third-order  error-sampled  system  with  zero-order-hold  and  tacho- 
meter feedback  shown  in  figure  2.   Suppose  that  there  is  a  requirement  that 
roots  for  the  sampled  transfer  function  have  zeta  values  greater  than  0.5  or 
lie  within  the  corresponding  spiral  contour  in  the  z-plane.  Allowable  s- 
plane  and  z-plane  root  areas  are  shaded  in  figures  3a  and  3b  respectively. 
Proper  adjustments  of  k  and  k  to  meet  this  requirement  are  available  from 
interpretation  of  the  parameter  plane  image  of  the  z-plane  contour  which 
bounds  the  allowable  roots.  Choose  oL    and  /3   as  k  and  kk  respective- 
ly; for  this  choice  the  characteristic  equation  of  the  system  shown  in  figure 
2  is: 

f(  od  ,  A  ,  *)  -  z3  +  (.058  oL  +  .310(3   -  1.974)z2  +  (32) 

(.16306  -  .122  j3  +  1.198)z  +  (.Q280<L  -  .188/3   -.223)  =  O 

The  coefficients  of  this  characteristic  equation  are  linear  in  OC    and  3   ; 

o£—  G    plane  curves  (mappings)  are  obtained  by  utilizing  the  parameter  plane 

formulas  for  the  linear  coefficient  case  developed  in  chapter  2.  These 

image  curves  and  the  images  of  the  two  points  at  which  the  designated  z- 

plane  contour  intersects  the  real  axis  are  shown  in  figure  4.   Shading  in 

this  figure  is  by  inspection.   (In  this  case  inspection  shading  can  be 
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achieved  as  follows.   Recognize  that  the  characteristic  equation  is  of  order 
three  in  z;    further  recognize  that  for  o£-  j3  -   0  the  roots  to  the  character- 
istic equation  are  identical  to  the  open  loop  poles  which  are  known  to  be 
within  the  contour.   Hence  all  three  roots  to  the  characteristic  equation  lie 
within  the  designated  contour  for  M-points  located  in  the  region  which  in- 
cludes the  parameter  plane  origin.   Since  this  is  the  maximum  number  of  roots 
which  can  be  enclosed  an  M-point  traveling  across  any  boundary  to  this  region 
implies  roots  leaving  the  chosen  z-plane  contour.  With  these  observations 
shading  can  be  effected  without  recourse  to  formulas).   Once  the  para- 
meter plane  curves  are  drawn  and  shaded  (as  in  figure  4)  the  compensation 
problem  reduces  to  locating  the  M-point  in  the  region  which  results  in  all 
roots  lying  within  the  specified  z-plane  contour.   The  preceding  discussion 
indicates  that  this  region  is  the  one  which  includes  the  parameter  plane 
origin.   To  verify  this  analysis  six  M-points  in  four  different  enclosed 
root  regions  have  been  indicated  in  figure  49   and  the  roots  to  the  character- 
istic equation  for  each  choice  have  been  computed.  Table  2  lists  these  roots 
and  indicates  how  many  fall  within  the  specified  z-plane  contour.  For  each 
M-point  theory  and  practice  agree. 

Let  it  be  emphasized  that  this  example  illustrates  parameter  plane  inter- 
pretation and  shading  procedures  rather  than  servo-system  design  in  the  para- 
meter plane.   For  the  latter  more  contours  (images)  are  drawn;  shading  may 
or  may  not  be  required  and  additional  information  to  guide  selection  of  an 
M-point  may  be  displayed. 

Mitrovic  Curve  Interpretation  for  Common  Z-Plane  Contours 

Mitrovic  methods  (defined  in  chapter  2)  are  special  cases  of  the  para- 
meter plane  which  warrant  attention  in  linear  (s -plane)  theory  because  the 
Cauchy  criterion  regarding  argument  change  for  roots  enclosed  by  s-plane 
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contours  is  directly  interpretable  as  encirclements  of  the  M-point.   For 
the  Mitrovic  plane  the  tedious  sub-division  of  the  parameter  plane  into  accept- 
able and  inacceptable  regions  can  be  bypassed  in  favor  of  a  simple  inspection 
which  determines  whether  the  Cauchy  criterion  for  enclosed  roots  is  or  is  not 
complied  with.   The  relation  between  constant  zeta  curves  on  the  Mitrovic 
plane  and  the  Cauchy  criterion  has  been  thoroughly  explored  in  the  literature. 
(6,  8).  A  compendium  of  these  relations  is  simply  this.  When  all  roots  of 
the  characteristic  equation  are  enclosed  by  the  s-plane  contour  being  mapped, 
argument  f(s)  where  f(s)  becomes  the  characteristic  equation  when  equated  to 
zero,  increases  monotonically.   In  turn  this  requires  that  the  M-point  be  en- 
circled in  a  predictable  manner  by  the  Mitrovic  plane  image  of  the  s-plane 
contour.   Were  constant  zeta  contours  meaningful  in  the  z-plane,  s-plane 
interpretation  would  be  applicable  to  them;  however,  interesting  z-plane 
contours  are  spirals  of  constant  damping  and  circles  of  constant  settling 
time.   It  can  be  shown  that  when  all  roots  of  the  characteristic  equation 
are  within  either  of  these  contours  argument  f(z)  increases  monotonically 
as  the  contour  is  traversed  in  a  positive  direction.   (Proof  can  be  obtained 
by  resolving  f (z)  into  its  factors  and  showing  that  each  has  a  monotonically 
increasing  argument  if  it  is  located  within  the  contour).  At  this  point 
complete  analogy  between  s-plane  and  z-plane  interpretation  of  Mitrovic 
curves  fails.   The  monotonic  increase  in  argument  still  sllows  satisfaction 
of  the  Cauchy  criterion  to  be  determined  from  the  manner  in  which  image  curves 
encircle  the  M-point  in  the  Mitrovic  (B  -B.)  plane.   However,  a  few  simple 
examples  suffice  to  show  that  this  interpretation  does  not  extend  to  general- 
ized Mitrovic  coefficients  for  z-plane  contours  of  constant  damping  and  con- 
stant settling  time* 
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Recall  that  the  characteristic  equation  and  Mitrovic  polynomial  are 
defined  as: 

f(z)  =  F(z)  +  a^  +  aQ 

f  (z)  =  F(z)  +  B  (z)z  +  B.(z)  =  0 
m  10 

Subtraction  of  zero  (the  Mitrovic  polynomial)  from  the  first  equation  yields: 

f(z)  =  (3l  -  B1(z))z  +  <a0  -  B0(z)) 

Consider  the  f(z)  locus  and  its  Mitrovic  (B--  B  )  curve  shown  in  figures  5a 
and  5b  respectively. 
Im[f(z)] 


Re[f(z)] 


Figure  5a  Figure  5b 

If  all  zeros  of  f(z)  are  enclosed  by  the  z-plane  contour  the  f(z)  locus  al- 
ways starts  on  the  positive  real  axis  and  argument  f(z)  increases  monotoni- 
cally  as  the  locus  f(z)  encircles  the  origin.  The  total  change  in  argument 
f(z)  is  IT      radians  times  the  number  of  roots  of  f(z)  enclosed  by  the  z- 
plane  contour.   Since  argument  f(z)  is  a  continuous  function  (assuming  no 
root  of  f(z)  lies  on  the  specified  z-plane  contour)  the  Cauchy  criterion  re- 
quires that 

argument  f(z.)  =kT     k=l,  2,  3,  ,  n-1 
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(that  is  f  (z.)  be  real)  for  n-1  values  of  z  on  the  specified  contour  where 
f(z)  is  a  polynomial  of  degree  n.   the  n-1  values  of  z  do  not  include  the 
two  points  at  which  the  specified  contour  intersects  the  real  axis.  At 
these  two  points  f(z)  must  also  be  real  but  the  interpretation  which  follows 
will  not  be  valid.   Take  particular  note  that  under  these  restrictions  the 
imaginary  part  of  z.  is  not  zero  at  these  n-1  points.   Hence 

f(z)  =  [ax    -  B^z)]  z.  +  [aQ  -  B^z.)] 

being  real  requires  that  [a.-B-Cz  )]z.  be  real.   This  occurs  only  if  a  - 

B,(z.)  =  0;  which  in  the  Mitrovic  plane  requires  that  the  B.-B,  curve  lie 
1   i  r       i  0   1 

on  the  vertical  line  through  the  M-point.  Furthermore  since  argument  f(z) 
is  monotonically  increasing  the  half  lines  a>B  and  a  <^Bn  must  be  cut 
alternately  with  the  half  line  an  <  B.  cut  first.  That  is: 

argument  f  (z)  =  2k IT         implies  a~  >B 

argument  f(z)  ■  (2k  +  1)7T    implies  aQ  <B 

The  half  line  an  <CB_  must  be  cut  first  because  argument  f(z)  is  monotonic 
and  must  equal  7T  >  a  point  on  the  half  line  an<B_,  before  it  can  equal 
2  IT  ,  the  point  at  which  it  crosses  the  half  line   an>  B~. 

The  points  at  which  the  specified  z-plane  contour  intersects  the  real 
axis  were  excluded  because  previous  arguments  did  not  apply  to  them.  Observe 

that  if  z.  is  real 

l 

f(z.)  =  [ax   -B^z.)]  z.  +[aQ  -B0(z.)]  i*  real 

and  there  is  no  requirement  that  [z  -B  (z  )]  =  0.   In  general  this  term  is 

not  equal  to  zero  for  real  z  and  this  fact  can  be  used  to  determine  the  manner 

in  which  the  Brt  -B,  curve  encircles  the  M-point.   If  z  is  the  value  at  which 
0   1  r 

the  specified  z-plane  contour  intersects  the  positive  real  axis  then: 
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B  C'z  )  2>  a        implies  counter  clockwise  encirclement 
B..(z  )  <^   a        implies  clockwise  encirclement 

This  follows  directly  from  the  requirement  that  the  line  a   <  B_  must  be 
cut  first.   The  value  B. (z  )  can  be  computed  by  the  relatively  simple  formula 

A/ 

Bl(Zr)  =  "2  KaKZr" 
Mitrovic  B  -  B  curve  interpretation  for  z-plane  contours  can  be  summarized. 

1.  If  all  roots  of  f(z)  are  located  within  a  specified  z-plane 
contour  (settling  circles  or  damping  spirals)  the  B  (z)  - 
B.  (z)  locus  encircles  the  M-point  in  the  Mitrovic  plane. 

2.  The  direction  of  encirclement  is  determined  by: 

B,(z  ) 2>  a,    implies  counter  clockwise  encirclement 
1  r     1 

B.  (z  ).<  a     implies  clockwise  encirclement 

3.  The  line  aft  -  B  =  0  must  be  intersected   (n-1)  times  plus 
one  additional  time  each  if  the  B_  -  B  curve  starts  or  ends 
on  this  line. 

In  linear  systems  the  next  attempt  was  to  extend  Mitrovic' s  method  to 
include  any  two  coefficients  of  the  characteristic  equation.  Desirable  regions 
in  the  generalized  plane  were  again  found  by  inspection  in  the  same  manner 
as  for  the  B  -  B.  case.  No  such  simple  interpretation  exists  for  z-trans- 
forms  when  generalized  Mitrovic  coordinates  are  used.  Since  zeta  varies  in 
all  z-plane  contours  of  interest  the  requirement  that  the  B.-B.  locus  en- 
circle the  M-point  (a. -a.)  is  no  longer  applicable,  (i.e.,  the  Cauchy  criter- 
ion can  no  longer  be  interpreted  in  the  B.-B.  plane  by  inspection).  A  few 
illustrations  will  suffice  to  illustrate  the  point. 

Illustration  I:     f (z)=  IV-B^z)!  z  +  Tag-B^z)-! 
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The  criterion  "f(z)  real"  is  met  when  z  =  j  (without  restrictions  on 
the  terms  in  brackets).   Thus  alternate  cutting  of  the  half  lines  B  "> 
a  and  B  <a  is  not  required  and  the  M-point  ne6d  not  be  encircled 
to  meet  the  requirement  that  argument  f (z)  be  monotonically  increasing. 
Illustration  II:    f(z)  =  [  {ak+1"Bk+1<z)J  z+  lvVz>3  ]  Z 
The  criterion  "f (z)  real"  is  satisfied  if  argument  z  equals  minus 
argument   [  V a,  .  -B,  1  (z)  ?  z  +  ^a, -B,(z)  j  J.  Again  the  argument 
requirement  on  f(z)  cannot  be  determined  by  simple  inspection  of  the 

B,  , ,  -B,  curve . 
k+1  k 

Although  it  is  probable  that  criteria  for  desired  Mitrovic  plane  regions  for 
generalized  coefficients  can  be  found,  they  will  be  complicated  and  it  is 
recommended  that  parameter  plane  interpretation  methods  be  used.  And*  with- 
out simple  interpretation  methods  the  Mitrovic  plane  loses  its  special 
attraction  over  the  parameter  plane  and  reverts  to  simply  a  special  case  of 
the  latter. 
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CHAPTER  IV 

Preceding  chapters  endeavored  to  present  and  interpret  the  character- 
istic equation  in  a  meaningful  manner  in  the  parameter  plane.   Parameter 
plane  presentations  need  not  be  restricted  to  the  characteristic  equation; 
in  a  more  general  consideration  the  parameter  plane  is  a  space  in  which  a 
system  is  described  in  terms  of  its  two  variable  parameters.   In  analogy 
to  characteristic  equation  presentation,  system  description  in  the  para- 
meter plane  is  the  reverse  of  descriptions  normally  seen.   Parameter  plane 
presentation  displays  suitable  values  of  oC   and  (3   to  meet  a  specified 
system  performance  criterion  whereas  conventional  procedures  first  fix  ol 
and  (3  and  display  system  performance  resulting  from  this  choice.   Signi- 
ficant advantages  acrue  when  system  criteria  are  described  in  the  parameter 
plane. 

1.  The  usual  "trial  and  error"  search  of  parameter  settings  to 
meet  specifications  is  eliminated.   This  advantage  is  espe- 
cially attractive  when  specifications  are  of  the  form  "equal 
to  or  less  (greater)  than"  and  as  such  are  not  tractable 
mathematically.   In  these  instances  the  locus  satisfying 
equality  usually  divides  the  parameter  plane  into  accept- 
able and  unacceptable  regions  where  parameter  settings  can 
be  chosen  by  inspection  to  satisfy  both  this  and  other 
criteria. 

2.  Simultaneous  consideration  of  all  criteria  and  character- 
istic equation  root  information  is  accomplished  by  inspec- 
tion.  Tacitly  assumed  is  that  each  criterion  plots  as  a 
curve  which  divides  the  parameter  plane  into  regions  suit- 
able or  unsuitable  for  that  criterion.   Parameter  settings 
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within  the  intersection  of  all  suitable  regions  simultaneous- 
ly meet  all  specifications. 

Liabilities  are  also  present  in  parameter  plane  presentations  although  they 

might  more  properly  be  classified  as  limitations. 

1.  System  descriptions  are  limited  to  description  in  two  para- 
meters (even  though  additional  parameters  may  be  available 
for  adjustment) . 

2.  Specified  criteria  may  be  difficult  or  impossible  to  deter- 
mine as  functions  of  adjustable  system  parameters.   A  maxi- 
mum overshoot  specification  typifies  such  criteria.   Theoreti- 
cally a  locus  of  o(^-j3     plane  points  which  have  the  specified 
maximum  overshoot  exists  but  in  practice  the  locus  is 
virtually  impossible  to  determine.   Time-integral  criteria 
and  magnitude  of  frequency  response  at  resonance  also  seem 

to  fall  within  this  class. 
These  limitations  not  withstanding,  two  very  important  performance  specifi- 
cations, steady  state  error  and  bandwidth,  can  be  displayed  on  the  para- 
meter plane.   They  complement  parameter  plane  information  already  avail- 
able from  characteristic  equation  roots  by  relating  system  performance  to 
the  zeros  of  the  closed  loop  transfer  function  as  well  as  to  the  poles. 
(The  formal  consideration  of  closed  loop  zeros  is  almost  universally  neg- 
lected by  other  transform  analysis  techniques) . 

Steady  State  Error  Curves  in  the  Parameter  Plane 

The  steady  state  error  or  inability  of  a  servo-system  to  follow  an 
aperiodic  input  disturbance  is  an  important  criterion  and  frequent  system 
specification.   For  linear  systems  or  linear  sampled  data  systems  steady 
state  error  is  usually  obtained  by  application  of  the  final  value  theorem 
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to  the  LaPlace  or  z-transform  of  the  error  function.   This  theorem  also 
allows  the  requisite  steady  state  error  curves  to  be  computed  for  para- 
meter plane  display.   Derivation  of  the  equation  governing  these  curves  in 
the  parameter  plane  will  be  for  the  z-transformation  with  a  subsequent 
paragraph  indicating  modifications  required  to  extend  applicability  to  the 
LaPlace  transformation  and  continuous  systems. 

Assume  that  the  open  loop  transfer  function  of  a  unity  feedback  system 
(or  a  system  which  has  been  reduced  to  unity  feedback  form)  is 


No(z) 


(z-1)-?  Do(z) 

where  D  (z)  contains  no  factors  of  (z-1)  nor  factors  in  common  with  N  (z) . 
o  o 

Further  assume  that  both  N  (z)  and  D  (z)  contain  no  terms  above  products 

o        o 

in  the  parameters  oC  and  8  .      These  assumptions  are  illustrated  in  figure 
6  and  equations  (33) . 


Mil 


V2> 


C(z) 


Figure  6 


Vz)  ■  Z  (*£+■  &  %  +  «fi  9J +  h\  )  z 


Dq(z)  =  ^  (c*  bK  4-  fi  CK   +  od/3c/K  4-  eK)  Z 


(33) 


K^O 

Typical  steady  state  error  specifications  require  aperiodic  inputs  of  steps, 
ramps,  accelerations,  whose  z- transforms  have  the  form 

-1 


R(z)  - 


where  R*(z)  contains  no  factor  of  (z-1)  or  (z-1) 


m  =  1  implies  step  inputs,  m  =  2  implies  ramp  inputs  and 
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m  =  3  implies  acceleration  inputs. 
The  z-transformed  error  signal  to  which  the  final  value  theorem  will  be 
applied  is  obtained  by  block  diagram  manipulation. 

«-jf  Do(z)R*(z) 


E(z)  =  (2-1) 


A/o(z)+  (2-0MDo(z) 


The  final  value  theorem  then  gives  the  steady  state  error  as 


)i 


E     =  II" 
s.s.    2.  — 


Do  ( z)  R*« 


(ZzJi  .  


(34) 


Consideration  of  equation  (34)  is  divided  into  three  cases. 

1.   jP+l-m^O  In  this  event  the  factor  of  (z-1)  or  (z-l)~ 

present  approaches  zero  or  infinity  as  the  limit  is  passed  and 
error  specifications  have  no  meaning  for  the  input  applied. 
Jf  +  1  -  m  =  0  and  X  =  0  This  denotes  a  type  zero  system  in 


2. 


which  the  steady  state  error  for  a  step  input  is 

DoO)  R*0) 


(35) 


3. 


X  +  1  -  m  =  0  and  X  £  0     This  denotes  type  I  systems  and  above 
in  which  a  factor  approaching  zero  as  a  limit  multiplies  D  (1) 


(36) 


in  the  denominator.   The  steady  state  error  expression  appropri- 
ate here  is 

E    .  p.  o)  mo 

Although  their  present  statement  does  not  so  indicate  equations  (35)  and 
(36)  are  the  relations  which  permit  steady  state  error  specifications  to 
be  drawn  on  the  parameter  plane.  Statements  adapted  to  parameter  plane 
use  can  be  obtained  by  evaluating  D  (1)  and  N  (1)  in  equations  (33)  and 
defining  a  more  concise  notation  for  the  resulting  sums.  That  is: 
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■.<»  ■  i  (* p* +■ p  *»  +■  "v*  % + hK ) 

MO 

Do(1)   =    2  (***  +  ?**  +  ^^  *  ^) 


Now  define: 

V 

■ti 

V 

V 

K=o 

V 

■S< 

G 


K=0 

Er.  "  J  £  (3?) 


K 


E 
s.s. 


With   these  definitions  N   (1)    can  be  expressed  as 
No(l)    =   C*Ep+    (5Eq+      ^EG+ER 
Do(l)   =    06  ER  +    /3  Ec  +      <yy8  ED  +  EE  (38) 

where  all  E*«  are  constants.   Useable  expressions  for  steady  state  error 
curves  on  the  oC- /3     plane  are  rendered  by  substitution  of  expression  (38) 
into  equations  (35)  and  (36).   For  substitution  into  (36) 

.  6*  CB +££+<*£  ED +  EV)R*P) 
(*E',**E,*VEft+Ei) 

In  this  expression  all  E's  are  constants  as  is  R*(l).   When  E     is  equated 

s  •  s  • 

to  the  steady  state  error  specification  the  equation  describes  the  locus  of 
all  c<-/8   plane  points  exactly  meeting  the  requirement.   The  locus  explicit- 
ly expressed  for  fi  as  a  function  of  oL   is 

A.       *  [£b  R*0)  -  EP  ltd  +  fee  R*U)-Eh  Ci.il  (39) 

*  Qeo  R*0)  -£*  F4S,]  +  [Ec  R"0>  -  E?  Es.  s] 

For  type  I  systems  and  above  equation  (39)  defines  the  locus  of  all  points 
(0Ct/3  )  having  a  steady  state  error  E     to  an  input  whose  z-transform 

I  S  •  S  s 
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is     R(z)   . 

For  type  zero  systems  an  analogous  formula  which  describes  steady  state 
error  loci  for  an  input  step  of  magnitude  R  can  be  obtained  by  initially 
substituting  expressions  (38)  into  equation  (35) .   The  final  result  is 

0  -  U  [gg  (R-Es-s.)  -EpEj.J  +  [f  »(*-£«.)  -  Eh  E*J 
°c[ED(ff-Es.s.)  -  Eo  £«.«]  4-  [Ec  (R-E&O  -  E<p  Es.3 
Observe  that  the  steady  state  error  loci  given  by  equations  (39)  and  (40) 
describe  hyperbolas  in  the  parameter  plane  for  the  oi/3   product  case. 

For  the  linear  case  E^  .!=  E„  =  0  and  the  steady  state  error  loci  become 

D1   G 

straight  lines  in  the  parameter  plane. 

Steady  State  Error  Loci  for  LaPlace  Transformed  Variables 

If  anything  steady  state  error  loci  are  easier  to  obtain  from  error 
signals  described  by  their  LaPlace  transform  as  opposed  to  signals  des- 
cribed by  their  z-transform.  As  for  the  sampled  data  case  assume  a  known 
open  loop  transfer  function  with  numerator  and  denominator  linear  in  o<L 
and  Q   but  not  necessarily  linear  with  respect  to  each  other.  That  is: 

N 

N0(s)  -2^PK°+/ff^+°^ar  +  /'^S  (41) 

do(.)  =  J (oc b°+-  ft c°K  4- oc/S d*  +  e°K )  sK 

where  D  (s)  does  not  contain  s  as  a  factor  and  the  complete  denominator 
o 

18  s^  D  (s) 

o 

Again  consider  only   inputs  of  the   step,    ramp,    acceleration  type. 

Input  =  Rt*"1  i^input]      =    fazllLg 

o 

Under  these  assumptions  and  definitions  the  LaPlace  transform  of  the  error 
signal  is: 

E(s)  -      S*Do(s)        m  (m-i)!  R 
/V0(s)  +  s*D0(s)        5Qs^ 


For  stable  systems  the  final  value  theorem  can  be  applied  to  exhibit  steady 
state  error  as: 


s.s. 


S-»o  L      J     i 


S-*o 


J+i-n      (Ai-i)/  R0o(s) 


A/0CS)-hS^00(S)_ 


(42) 


On  passing  the  limit  results  can  again  be  divided  into  three  cases. 

1.  X+   1-m  ^  0  In  this  case  steady  state  error  is  either  zero  or 
infinite  and  no  finite  adjustment  of  the  parameters  can  alter 
this  result. 

2.  X  +  l"m  =  Of   Jir  0  These  requirements  describe  a  system  of  type 

I  or  higher.   The  parameter  plane  locus  of  points  having  a  steady 

state  error  of  E      as  obtained  from  equations  (41)  and  (42) 
s.s. 

is: 

_    o4;  (ai-i)/  r  - p:es2+& (w-q//?  -  a;  es,j       (m) 

3.  ,\+  1-m  =  0,  Jc  ■  0  This  requirement  is  indicative  of  type  zero 

systems  which  are  characterized  by  having  steady  state  errors  to 

step  inputs.   The  parameter  plane  loci  for  a  steady  state  error 

E     to  an  input  step  of  magnitude  R  are  governed  by  the  equa- 
s  •  s  • 


tion: 


/8 


<*  [g.  ( R-  Et.Q  -ft-gJ  +  jk  f  R-e»)  -hi  eJ 


<*  [«crt-E„)  -?;£•«]  + 


(44) 


_c:(r-e„)-?;£sj. 

Observe  that  application  of  the  final  value  theorem  makes  only  those  terms 
in  the  numerator  and  denominator  not  multiplied  by  s  (after  s  factors  in 
numerator  and  denominator  are  removed)  significant  in  determination  of 
steady  state  error.   If  parameters  oC  and  /6  do  not  appear  in  these  terms 
steady  state  error  is  not  influenced  by  their  adjustment  regardless  of  where 
they  appear  elsewhere.   Furthermore,  only  the  manner  in  which  oC  and  |S 
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appear  in  these  terms  affects  the  shape  of  the  steady  state  error  curves. 

If  o6fi  products  appear  in  these  terms  the  steady  state  error  loci  will 

be  hyperbolas;  if  oC   and  /G  appear  linearly  the  loci  will  be  straight  lines. 

Bandwidth  Curves  in  the  Parameter  Plane  * 

The  bandwidth  is  significant  because  it  indicates  rise  time  or  speed 
of  response,  it  measures  in  part  the  ability  of  the  system  to  reproduce 
the  input  signal,  and  it  approximately  describes  the  filtering  character* 
istics  of  the  system.1 

Bandwidth  specifications  can  be  presented  on  the  parameter  plane  with  no 
restrictions  above  those  inherent  to  explicit  solution  of  OC    and  /3 
When  so  displayed,  bandwidth  can  be  considered  in  conjunction  with  other 
specifications  presentable  on  the  parameter  plane  and  satisfaction  of  all  is 
virtually  by  inspection,  if  such  satisfaction  is  possible.  Before  embarking 
on  the  derivation  of  formulas  which  display  bandwidth  criteria  on  the  para- 
meter plane  some  preliminaries  deserve  attention.  Among  these  is  agreement 
on  a  definition  for  bandwidth. 

Excite  a  servo-system  by  a  sinusoidal  test  signal  and  form  the  ratio  of 
peak  fundamental  component  output  to  peak  input.   (If  the  system  is 
linear  the  output  consists  wholly  of  a  fundamental  component).   The 
lowest  frequency  at  which  this  ratio  falls  below  a  designated  value  A 
is  defined  as  the  bandwidth  of  the  servo-system. 
The  most  widely  accepted  value  for  A  is  0.707  (-3  db)  although  unity  and 
other  values  have  been  used.  A  readily  accessible  expression  from  which  the 
ratio  of  peak  fundamental  output  to  peak  input  can  be  computed  is  the  closed 
loop  transfer  function.   If  Fourier  transforms  are  employed  the  magnitude 
of  the  closed  loop  transfer  function  represents  this  ratio.   The  identical 
ratio  is  obtainable  from  the  magnitude  of  the  closed  loop  transfer  function 

Truxal,  J.  G.  Control  Systems  Synthesis.  McGraw-Hill,  1955:  77. 
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expressed  in  terms  of  the  complex  LaPlace  variable  s  if  the  substitution 
s  =  jw  is  effected.   For  sampled  signals  a  similar  ratio  is  obtained  by 
examining  the  closed  loop  z-transfer  function  for  values  of  z  on  the  unit 
circle. 

Parameter  plane  methods  are  not  primarily  concerned  with  a  determina- 
tion of  this  ratio;  rather  the  values  of  system  variables  (locus  of  oC-^ 
plane  points)  which  meet  bandwidth  specifications  (cause  the  output-input 
ratio  to  be  of  proper  magnitude  at  the  specified  frequency)  are  sought. 
Functionally  the  parameter  plane  bandwidth  problem  is  expressed: 

(45) 


A^ 


R(uJew) 


where  by  adjustment  of  system  parameters  °C   and  /3  the  magnitude  of  the 
closed  loop  transfer  function  is  forced  to  have  a  magnitude  A  at  the  band- 
width frequency  w    .  Equation  (45)  while  sufficient  to  define  bandwidth 
loci  on  the  parameter  plane  is  not  conveniently  utilized  for  their  construc- 
tion. The  remainder  of  the  section  is  devoted  to  remedy  of  this  problem. 

In  line  with  parameter  plane  assumptions  the  closed  loop  system  trans- 
fer function  is  required  to  have  the  form 

Kz'  _  JL=° (46) 

R(z)~ 


2  (^K  +  zScK+^/H+eJz 


K=0 

In  the  course  of  the  derivation  it  proves  convenient  to  separate  both  the 
numerator  and  the  denominator  into  their  real  and  imaginary  components. 
The  Chebyshev  functions  introduced  in  equations  (5)  through  (9)  provide  a 
satisfactory  vehicle  for  this  purpose.  For  reference,  z  defined  in  terms 
of  the  Chebyshev  functions  is 


To  promote  brevity  the  real  and  imaginary  parts  of  each  sum  are  defined  as: 

v  5.  •  u<*  >  -  2  *« ^  f-^  u*(0 

V 

The  remaining  twelve  quantities  in  this  set  are  similarly  defined.  It  is 

noteworthy  that  by  this  definition  all  quantities  R  and  I  are  real.  On 

replacing  the  sums  in  equation  (46)  by  their  defined  equivalents  in  equa- 
tion (47)  one  obtains: 

c(i)  m  oc  (*P  *- ix£±££a± i ?<,)+*/%  +3^)»(*m-aii.)     (48) 

R(Z)         u.  (R8-h  Jl6)+  /S(KC  +  a^  +  ^^JloH^r*-  Jlr) 
Bandwidth  specification  that  the  ratio  of  sampled  output  to  sampled  input 

have  some  magnitude  A  at  a  bandwidth  frequency  of  w,    equates  to  the  re- 

quirement  that  C(z,    )  over  R(z    )  have  magnitude  A  where  z     is  the 

O » W  •  D  •  W  •  D  •  W# 

value  of  z  on  the  unit  circle  which  images  the  bandwidth  frequency  s  ■  jw, 

D  ■  W  ■ 

in  the  s-plane.  When  these  bandwidth  conditions  are  introduced  the  quanti- 
ties R  and  I  become  constants  and  equation  (48)  describes  the  locus  of  all 
points  (oC  f  Q    )  which  satisfy  system  bandwidth  specifications.  That  is: 


A  = 


C(Zi.^) 


R(2k.O 

(«  f?B+  /QRC  ^RDf  RE)%  (*I8+*It  +  ^ID+Ie)2 
where  A  is  the  magnitude  of  attenuation  specified  at  the  bandwidth  frequency. 
All  R's  and  I's  are  constants. 

Performing  indicated  multiplications  allows  this  equation  to  be  render- 
ed as: 
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A2=  ojgfffi-^)  j-  e^/a^-Hpi^)  +  g*lfiM,+iQjfc)  +        49) 


Equation  (49)  describes  the  locus  of  parameter  plane  points  which  satisfy 
a  given  bandwidth  specification  demanding  an  attenuation  A  at  bandwidth 
frequency  w,    .   Its  appearance  is  forbidding  but  after  recognition  of 
certain  features  the  equation  itself  is  not  formidable.  Foremost  is  that 
c<  and  f3   individually  appear  to  no  power  higher  than  two.   Thus  to  find 
a  point  on  this  bandwidth  locus  for  a  given  alpha  value  involves  no  more 
than  solution  for  y6  by  the  quadratic  formula.   (Computer  plotting  of  such 
a  curve  is  trivial) .   It  is  also  significant  that  equation  (49)  describes 
a  conic  section  whenever  the  numerator  and  denominator  coefficients  are 
linear  in  o^  and  f$   in  which  case  R  =  I  =  IL=  I  =0. 

To  demonstrate  the  first  of  these  observations  explicit  solution  for  13 
will  be  effected  under  the  assumption  that  oL   is  known  (assigned).   For 
convenience  first  define: 

Xx(  d.  )  -  0<2  <R2  +  !*)  +  2  *  (RQRG  ♦  IQIG)  +  (R2  +  I2) 

X2<<*  )  =  c*2  <l£  +  I*)  +  2  OC«ch  *  lch)   +  (^  +  !2) 

Y^OC)   =   2oC2(RpRG  +  IpIG)  +  2«V     +  R^  +  IT     +  IGIH)  + 


2<RQRH  +  1QIH> 


(50) 


55 


Y2(od  )  =  2  ot2^  +  IBID)  +  2  oC^  +  %RE  +  IBIC  +  IdIe)  + 

2<Ve  +  ICIE) 
Wx(oC  )   =   oc  (^  +  IP)  +   2  °^<RpRH  +  IPIH)  +   (RH  +  4 


W2(o<l  )   =  OCZ(*l  +  Ig)  +  2  <X  (R^  +  ibIe)  +  R^  +  I*) 


Explicit  solution  of  /?  from  equation  (49)  utilizing  the  notation  defined 
in  (50)  is:  

P°  ?(^Xe-X,)  <51) 

Equation  (51)  permits  determination  of  the  jg  coordinate  of  parameter  plane 
points  meeting  a  bandwidth  specification  given  the  value  of  the  oC  coordin- 
ate. A  negative  value  under  the  radical  sign  infers  that  no  point  /3 

will  satisfy  the  bandwidth  specification  for  the  value  of  oC     assigned. 

w 
When  bandwidth  is  specified  at  a  frequency  of  w  ■  s/„  the  bandwidth  loci 

in  the  parameter  plane  degenerate  from  conic  sections  into  straight  lines 
(for  the  linear  coefficient  case).  Whether  interpretation  of  these  degener- 
ate curves  is  consistent  with  interpretation  of  bandwidth  curves  specified 

at  other  frequencies  has  not  been  established;  hence,  the  reader  is  caution - 

w 
ed  to  avoid  specifying  bandwidth  at  the  frequency  of  s/9« 

Bandwidth  Curves  for  S -Plane  Transfer  Functions 

The  preceding  bandwidth  theory  is  applicable  to  closed  loop  transfer 
functions  described  in  terms  of  the  complex  variable  s.   In  fact  only  the 
definitions  of  the  R's  and  I's  need  be  altered  to  make  z-plane  bandwidth 
equations  valid  for  s-plane  use.   For  the  sake  of  completeness  the  deriva- 
tion below  duplicates  some  of  the  work  in  the  preceding  section.   Assume  a 
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closed  loop  transfer  function: 


CCS)  . 


Cft)   £  (*  Pk  »  /»»»  -  °^ ?«  + hK)LOt   |T«  (- CK  3 F?UK(- C )] 

Since  bandwidth  is  always  specified  at  some  real  frequency  only  zeta  values 
of  zero  are  relevant.  For   =  0  the  T.  (0)  and  U.  (0)  are  expressible  as: 


T.  (0)  »  cos  JUL  n  (0)  =  sin 


KIT 


kN  g  k  g 

This  observation  allows  definitions  of  the  R*s  and  I's  in  equations  (47)  to 
be  altered  for  s-plane  use. 

W  =Z/"^  9**  *  W  =^_/)  9«+'Wn   (52) 

K*0  *=° 

Omitted  R's  and  I's  are  analogously  defined. 

S-plane  bandwidth  specifications  require  the  magnitude  of  output  to 
input  to  have  suffered  an  attenuation  A  at  the  bandwidth  frequency  w 
Functionally: 

|C(MJ| 


A  - 


R(dUL)\ 


Evaluate  equations  (52)  at  w  =  w,    and  substitute  the  result  into  the 
n  n    b.w. 


equation  above  to  obtain: 


A  - 


From  this  point  all  equations  (49)  and  subsequent  for  z-plane  bandwidth 
determination  apply  identically  for  s-plane  bandwidth  determination. 
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(53) 


, 


Periodicity  of  Bandwidth  Problem 

A  periodic  frequency  spectrum  is  characteristic  of  any  signal  described 
by  z-transforms.   Frequently  this  periodicity  characteristic  is  not  intrin- 
sic to  the  signal  being  described  but  is  artifically  introduced  by  a  fictit- 
ious sampler  to  facilitate  mathematical  analysis.   The  commonplace  situa- 
tion for  sampled  data  systems  is  illustrated  in  figure  (7). 


real 
sampler 

c(t) 

r(t)/~ 

system 

!   --'   I  c*(t) 

\ 

) 

1 

i 
i 
i 

__J 

\ 

fictitious 
sampler 

Figure  7 
This  frequency  spectrum  periodicity  present  in  z-transfer  functions  radi- 
cally alters  frequency  response  interpretation  and  invalidates  or  restricts 
prevalent  frequency  response  criteria  when  applied  to  pulsed  transfer  func- 
tions.  Among  criteria  so  altered  is  bandwidth;  the  3  db  attenuation  us- 
ually specified  for  servo-system  bandwidth  may  never  occur  when  testing 
the  pulsed  transfer  function  yet  the  system  itself  may  be  perfectly  satis- 
factory.  This  dilemma  can  profitably  be  bypassed  by  return  to  system  des- 
cription in  the  s-plane  even  though  such  description  involves  transcenden- 
tal functions.   Here  bandwidth  still  has  meaning  and  attention  need  not  be 
restricted  to  sampled  signals,  a  significant  advantage  when  examining  the 
system  output. 

Sampling  makes  the  previous  definition  of  bandwidth  somewhat  inade- 
quate (as  will  subsequently  be  shown)  and  agreement  on  a  suitable  defini- 
tion of  bandwidth  is  incumbent.  Consider  the  error  sampled  system  in 
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figure  (8). 
R(s) 


T 


C(s) 


Figure  8 
Denote  by  superscript  X*(s)  the  LaPlace  transform  of  a  sampled  function 
whose  continuous  transform  is  X(s).   Reasonable  alternatives  for  definition 
of  bandwidth  are: 


lc*(s)| 
U?*(s)| 


C(S)[ 

R(s)     A 


A  = 


|C(s) 


R*<fs)| 


s=ju; 


bw 


The  first  of  these  alternatives  is  identical  to  results  obtained  using  z- 
transfer  functions  and  its  inadequacy  has  already  been  discussed.   The 
second  of  these  alternatives  allows  a  different  bandwidth  for  inputs  R.f s) 
and  R~(s)  having  the  same  value  at  the  sampling  instants.   For  the  error 
sampled  system  shown  the  output  is  unaffected  by  the  value  of  the  input 
at  other  than  the  sampling  instants  and  one  suspects  that  any  definition 
of  bandwidth  should  also  have  this  property.   Therefore  bandwidth  appears 
most  suitably  defined  when  stated  as: 

Even  this  statement  needs  minor  modification  to  be  intuitively  satisfying. 
The  example  in  the  following  section  will  demonstrate  a  more  satisfactory 
and  inherently  pleasing  definition  for  bandwidth  is: 


A  -  — 


|CCS)| 


T     l**CS)|   w 


(54) 
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It  is  acknowledged  that  this  definition  for  bandwidth  is  subject  to  debate 
when  a  system  is  not  input  or  error  sampled. 

Block  diagram  manipulation  of  the  system  in  figure  (8)  shows  the  in- 
put to  be  related  by 

This  relation  allows  expression  of  system  bandwidth  as: 
G(S) 


T 


/4-  GH*(S) 


(55) 


A  proper  approach  toward  rendering  this  bandwidth  expression  fit  for  para- 
meter plane  use  might  be  to  examine  its  individual  terms.   G(s)  is  the 
system's  LaPlace  transfer  function  which,  with  the  possible  exception  of 

a  hold  element,  is  presumed  linear.   As  such  G(s)  must  be  a  ratio  of  two 

sT 
polynomials  in  "s"  and,  if  a  hold  element  is  included,  e   .   GH*(s)  is  the 

LaPlace  transform  of  a  sampled  function  and  it  can  be  represented  as  an  in- 
finite series. 


x*(s)  =^T  X(nT)€ 


zv=o 
where  x(t)  is  the  impulse  response  of  component  X(s).   Whenever  G(s)  and 
H(s)  are  linear  except  for  hold  circuits,  as  is  assumed,  the  series  for 

GH*(s)  has  a  closed  form  representation  which  is  a  rational  function  of 

sT 
e   .   Significantly  this  closed  form  is  readily  obtainable  from  z-trans- 

form  tables. 


GH*(s)  =  GH(z)      ST  (56) 

'z  =  e 

Also  noteworthy  is  that  specification  of  bandwidth  at  a  frequency  of  s  = 
jw     ■  jw  /2  forces  GH*(s)  to  be  real.   (i.e.   Specify  an  input-output 

0|W|  S 
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attenuation  of  magnitude  A  at  w  /2  in  preference  to  specification  of  an 

s 

attenuation  A'  at  some  other  frequency).   Such  specification  reduces  the 

problem  of  separating  the  real  and  imaginary  parts  of  the  function 

G(S) 

/+GH*(S) 
preparatory  to  obtaining  its  magnitude,  to  a  problem  of  separating  the  real 

and  imaginary  parts  of  G(jw)  which  is  a  rational  function  of  w. 

Still  requisite  is  an  equation  which  allows  ready  determination  of  para- 
meter plane  points  satisfying  the  bandwidth  specification  of  equation  (55). 
Assume  that  the  linear  or  product  appearance  of  oC  and  B    in  the  closed  loop 
transfer  function  continues  when  bandwidth  is  as  specified  in  equation  (55). 
To  elaborate,  grant  that  systems  having  closed  loop  transfer  functions 
shown  in  equation  (46)  will  retain  this  &ame  type  oC- R      dependence  under 

equation   (55) . 

CCS)     _    «.Vs)+/??Q(S)  +  o^f6 ($)  +  £($) 
R*(s)        <*  f8(s)  +•  fi  £CS)  +  ocfi  fD  (s)  +  iE($) 

Evaluate  the  functions  f  at  the  specified  bandwidth  frequency  s  =  jw 
and  define: 

Rp  =  Real  Part   [fp(jWB.W.^ 

Ip  -  Imaginary  Part  [fp(JwR  w  )]  (57) 

R,-  Real  Part   t^O^)] 

Ifl  -  Imaginary  Part   [fB<JwB#w<>)l 

All  R's  and  I's  omitted  above  are  similarly  defined.  Also  define: 

A'  ■  A*T  (57a) 

Parameter  plane  bandwidth  loci  for  bandwidth  as  defined  in  (55)  can  be  ob- 
tained using  definitions  (57)  and  equations  (49)  and  subsequent  which  were 
originally  developed  for  z-plane  operations. 
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c(w) 


R*(ou;) 


Illustration  of  Auxiliary  Parameter  Plane  Curves 

The  generalization  that  examples  always  serve  to  clarify  theory  is 
presumptuous  especially  in  this  instance  when  bandwidth  and  steady  state 
error  loci  are  already  known  to  be  conic  sections  and  straight  lines  respec- 
tively. Nevertheless  an  example  is  presented  and  surprisingly  the  general- 
ization holds.   In  fact  a  great  deal  can  be  learned  from  this  example  if 
digression  to  include  several  closed  loop  frequency  response  curves  is 
allowed.  Consider  the  third  order  sampled  data  system  illustrated  in 
figure  (2)  and  direct  attention  to  its  closed  loop  frequency  curves 
for  representative  parameter  settings  of 

a)  k  =  0.5,       k  -  1.0  figure  (8a) 

b)  k  =  0.1        k  -  1.0  figure  (8b) 

t 

These  curves  are  drawn  in  figure  (8) .   Immediately  obvious  and  disturbing 
is  the  fact  that  the  closed  loop  d.c.  (ws0)  response  of  this  feedback 
system  has  a  magnitude  differing  from  unity.   One  intuitively  expects  event- 
ual equality  of  output  and  input  when  a  type  I  system  is  disturbed  by  a  d.c. 
signal;  an  expectation  which  equates  to  a  frequency  response  of  unit  magni- 
tude for  w=0.  The  cause  for  fallacy  of  intuition  and  a  definition  modifica- 
tion which  results  in  more  pleasing  (familiar)  frequency  response  magnitude 
curves  are  obtainable  from  examination  of  the  relative  magnitudes  of  R(s) 
and  R*(s).  To  expand,  the  magnitude  of  the  frequency  response  as  defined 
for  continuous  systems  and  also  the  definition  on  which  intuition  is  based 
is 

M  = 


C(4w) 


R(dw) 
whereas  a  wiser  choice  for  sampled  data  systems  seems  to  be 

C(:M 


M'  = 


R*(iej) 
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FIGURE  8a 
FREQUENCY  RESPONSE  CURVES 


FIGURE  8b 
FREQUENCY  RESPONSE  CURVES 
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Now  suppose  magnitude  R(jw)  does  not  equal  magnitude  R*(jw).  Then  correct 
intuitive  reasoning  precludes  the  same  result  for  sampled  data  systems. 
However  a  scale  factor  may  be  available  to  force  sampled  data  frequency 
response  curves  to  r6semble  their  counterparts  for  continuous  systems  more 
closely.  To  investigate  these  conjectures  consider  the  general  harmonic  in- 


put  e" 


and  form  the  ratio 


R*(S) 


i  W+)}  =  RCS)  =  ^ 


sT 


I 


jUrVa  =  R*(s)  -  e.T*e^T    ~   ,,e-(s-3uor 


PCs) 


-e 


-(5  -3U)0)T 


The  magnitude  of  this  ratio  is  of  interest  at  the  frequency  of  the  forcing 


function  only. 


R*huJo) 


/-e 


-is-i^T 


S-<lw0 


=  T 


Thus  an  impulse  sampler  attenuates  a  continuous  harmonic  input  signal  by  a 
factor  of  "T",  Such  attenuation  should  be  reflected  in  the  frequency  re- 
sponse curves;  a  fact  witnessed  by  the  example  curves  (solid i lines)  of 
figures  (8a)  and  (8b)  where  T  =  0.5.  Since  attenuation  factor  "T"  is  in- 
dependent of  forcing  frequency  w  ,  dividing  the  input-output  ratio  by  "T" 

o 

would  bring  the  zero  frequency  response  to  unity  and  the  sampled  data 
frequency  response  in  general  would  closely  correspond  to  that  for  continu- 
ous systems.   These  desirable  ends  supply  the  motive  for  the  bandwidth 
definition  of  equation  (55).  Observe  that  the  frequency  response  curves 
as  defined  for  z-transfer  functions 


A  = 


C(z) 


R(2) 


z  =  e 


no0r 


c*(s) 


/?*(s) 
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[shown  dashed  in  figures  (8a)  and  (8b)]  have  this  attenuation  factor  "T"  pre- 
sent in  both  numerator  and  denominator.  Hence  these  curves  are  conceptual- 
ly satisfying  without  recourse  to  scaling.  The  curves  of  figures  (9)  are 
like  those  of  figures  (8)  except  that  s-plane  frequency  response  definitions 

are  scaled  to  have  unity  value  for  inputs  of  w  =0  and  Y-axis  coordinates 

o 

are  in  decibels.  The  parameter  settings  for  two  of  the  three  curves  shown 

a)  k  =  0.44,      kfc  -  0.177         figure  (9a) 

b)  k  =  0.39,      kfc  -  1.643         figure  (9b) 

lie  on  the  same  parameter  plane  bandwidth  curve  in  figure  (10) ;  consequent- 
ly each  frequency  curve  should  exhibit  a  magnitude  of  0.3  (-10.5  db)  at 
frequency  w  =  -5-  .  Although  the  frequency  response  curves  of  figures  (9) 
apply  to  but  three  specific  examples  the  reader  should  recognize  the  extreme 
low  pass  properties  of  the  hold  element  and  the  close  correspondence  between 


T 


CCS) 


Rr(s) 


and 


C*CS) 


R*CS) 


for  w  <Z   -=£■   which  these  curves  illustrate.  These  observations  are  general- 
ly valid  although  the  frequency  at  which  the  generalizations  cease  tc  apply 
is  dependent  on  the  particular  system  and  the  sampling  rate. 

These  digressions  are  relevant  to  parameter  plane  study  to  the  extent 
that  they  allow  meaningful  specification  and  interpretation  of  bandwidth 
loci  thereon.   Several  such  loci  are  shown  (solid  curves)  in  figure  (10). 
For  the  specified  frequency  each  curve  divides  the  parameter  plane  into 
regions  having  attenuation  greater  than  or  less  than  the  attenuation  evidenced 
by  the  parameter  settings  on  the  curve.  If  the  curve  itself  represents 
points  just  meeting  a  bandwidth  requirement  then  parameter  settings  corres- 
ponding to  points  in  one  of  these  regions  will  all  exceed  the  bandwidth 
specification.  Moreover  no  intersection  of  bandwidth  loci  occurs  for  various 
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FIGURE  9  a 

FREQUENCY  RESPONSE  CURVES 
FOR  A  SAMPLED  DATA  SYSTEM 


ndtm  9b 


FREQUENCY  RESPONSE  CURVES 
FOR  A  SAMPLED  DATA  SYSTEM 
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FIGURE  9c 

FREQUENCY  RESPONSE  CURVES 
FOR  A  SAMPLED  DATA  SYSTEM 


The  bandwidth  and  steady  state  error  loci  displayed  in  figure  10 
are  for  the  third  order  sampled  system  shown  in  figure  2.  Equation 
55  governing  these  bandwidth  loci  is: 

<*  (l-e-sT) 


A=JL 


s^(s+l)(s+2) 


1  +  (.058*  ».3i/3  )e2sT  +  (.16*  -.12/?  )e3T  -H.026*  -,19/?  ) 
e3sT  -l.97e2sT4.1.20esT  -0.22 
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magnitude  specifications  at  a  given  frequency  although  curves  representing 
the  same  magnitude  attenuation  specified  at  various  frequencies  can  conceiv- 
ably have  common  points. 

The  dashed  straight  line  curves  in  figure  (10)  represent  steady  state 
error  loci.  Perhaps  their  most  significant  feature  is  the  vertical  slope 
from  which  the  reader  infers  that  steady  state  error  is  independent  of  the 
parameter   /S  ■  kk  . 
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CHAPTER  V 

Concepts  of  Dominant  Mode  Design 

The  preponderance  of  system  criteria  in  the  s-plane  are  based  on 
dominant  mode  analysis  and  design.   Such  concepts  as  damping  factor,  reso- 
nant frequency,  bandwidth  and  maximum  overshoot  are  meaningful  only  within 
the  context  of  a  dominant  mode  system.   (i.e.  What  meaning  has  bandwidth 
if  an  attenuation  of  -3  db  occurs  at  three  points  on  the  frequency  res- 
ponce  curve?)  Rudiments  of  dominant  mode  design  reveal  the  direction  any 
extensions  of  this  method  must  take. 

1.  An  ideal  model  must  be  defined.   (For  s-plane  use,  the  no- 
zero,  two  complex  pole  model  is  universally  accepted  as 
standard) . 

2.  The  actual  system  must  be  compensated  to  behave  like  the 
ideal  model.   (In  the  s-plane  the  actual  system  approximates 
the  ideal  second  order  model  if  a  pair  of  complex  poles 
dominate  its  behavior). 

Dominant  mode  design  in  the  parameter  plane  follows  the  indicated  path.  As 
an  initial  step  a  z-plane  "ideal  model"  is  defined  and  system  behavior  (rise 
time,  maximum  overshoot)  is  related  to  model  parameters.   The  following 
section  discusses  conditions  requisite  on  other  systems  to  insure  dominant 
mode  (or  model-like)  behavior.   Finally  the  parameter  plane  is  introduced 
as  a  vehicle  for  designing  or  compensating  the  actual  system  to  meet  the 
requirements  for  guaranteed  dominant  mode  performance. 

The  Dominant  Mode  Sampled  Data  System  Model  and  Its  Specification 

Choice  of  a  model  is  critical  because  dominant  mode  design  can  never  in- 
tentionally produce  systems  whose  properties  exceed  those  of  the  model.   (To 

illustrate  consider  an  s-plane  design  in  which  a  fourth  order  system  is 
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compensated  to  have  two  dominant  complex  poles.   It  may  be  that  the  fourth 
order  system  when  compensated  outperforms  the  second  order  model  but  such 
performance  is  coincidental  rather  than  designed).   An  elaborate  model  might 
be  justified  were  it  not  for  two  additional  requirements. 

1.  Common  specifications  must  be  readily  translatable  as  model 
parameters  (root  locations).   By  way  of  example,  a  maximum 
overshoot  specification  can  be  translated  into  a  zeta  re- 
striction on  the  roots  of  the  second  order  s-plane  model. 

2.  The  model's  pole-zero  configuration  must  have  simplicity 
sufficient  to  enable  the  designer  to  recognize  when  model 
like  behavior  (dominance)  is  achieved.   Suppose,  for  an 
example  by  contrast,  it  is  possible  to  determine  a  fourth 
order  ideal  model  system.   A  designer  now  compensates  an 
existing  system  to  approach  this  model;  but  approach  it 
where?   Should  he  design  to  have  two  poles  in  close  prox- 
imity and  two  further  removed  or  should  all  poles  be  moder- 
ately removed  from  those  of  the  model  or  is  some  other  ap- 
proximation to  be  preferred? 

These  considerations  indicate  an  acceptable  model  to  be  one  whose  closed 
loop  z-transfer  function  is 

C(z)  a     k(z--z,)  (58) 

where  z  is  a  real  zero 
1 

P^  an(*  p".  are  conjugate  complex  poles 

k  - .(/  -Pi)(  l~n)  This  choice  of  "k"  guarantees  a  type  I  system 

O-Z,) 
with  zero  steady  state  error  to  unit  step  inputs. 
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The  utility  of  this  model  is  contingent  upon  the  facility  with  which 
system  specifications  are  translatable  into  acceptable  pole-zero  configura- 
tions.  To  amplify,  unless  acceptable  model  pole-zero  configurations  can 
be  simply  determined  design  might  profitably  proceed  in  some  direction  other 
than  in  forcing  the  system  to  resemble  a  model.   The  following  paragraphs 
develop  the  required  relations  between  maximum  overshoot,  rise  time  and 
model  pole -zero  locations.   This  development  largely  duplicates  the  work 
on  pages  170  through  179  in  Kuo  and  is  included  to  accentuate  approximations 
inherent  in  the  model  itself  [2].   These  approximations  are  over  and  above 
those  made  when  compensating  higher  order  systems  to  a  dominant  mode  and  are 
not  required  for  continuous  systems. 

The  rise  time  and  maximum  overshoot  of  the  system  described  by  equation 
(58)  to  a  unit  step  input  are  desired.   This  requires  expression  of  c*(t). 

C(2)-(2-f?)(z-p,)  'T=T 

By  the  inversion  integral: 

c(nT)=-L-f KfZ-Z,)Z    Z"-'  J, 

r 

rT 
where  P    is  a  circle  of  radius  e  '  centered  at  the  origin  and  enclosing  all 

poles  of  C(z).   Applying  the  residue  theorem  to  this  integral  one  obtains 

0-PO-1B       (?-')(/»-?)       (f-D(f-n) 

which  under  definition  of 

1$ 


becomes 

c(nT),/-£k|g:f,K;j   s,N0r+M§) 
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or 


cfrn-<  +  8*-|gi*jffjjj      co,^  +  n$ 


where       V    -  arg^p^-z  )    -  arg(p  -1) 

-Gy  -  arg(p1-*1)  -  arg(Pl-l)  -  ^~ 
With  specification  of  the  closed  loop  poles  and  zero  the  above  expression 
is  sufficient  to  determine  the  output  c*(t)  at  the  sampling  instants.   How- 
ever this  expression  has  two  serious  defects. 

1.  Th*e  output  is  specified  only  at  the  sampling  instants  whereas 
peak  overshoot  in  the  continuous  output  is  likely  to  occur 
during  the  intersampling  period. 

2.  The  output  c*(t)  is  specified  as  a  function  of  three  para- 
meters, a)  real  part  of  complex  poles,  b)  imaginary  part  of 
complex  poles  and  c)  location  of  real  zero;  consequently  its 
universal  solution  is  difficult  to  display.   That  is  to  say  - 
system  criteria  which  are  functions  of  three  variables  are 
not  conveniently  represented  in  two  dimensions  or  less. 

The  next  step  is  circumvention  of  these  difficulties.   The  first  problem 

can  be  attacked  in  several  ways  among  which  are  employment  of  modified  z- 

trans forms  and  approximation  of  the  continuous  output  c(t)  by  a  continuous 

function  c  (t)  which  agrees  with  c(t)  at  the  sampling  instants.   The  complex- 
a 

ity  of  the  modified  z-transform  and  the  fact  that  zeros  of  C(z,m)  are  func- 
tions of  m  preclude  the  use  of  this  avenue.   Approximation  of  the  continu- 
ous output  with  a  continuous  function  which  agrees  at  the  sampling  instants 
offers  many  advantages.   Among  them,  the  peak  overshoot  of  this  approximating 
function  can  be  determined  and  if,  as  is  hoped,  the  approximation  is  good 
then  this  peak  should  be  close  to  the  true  peak  overshoot.   The  second  pro- 
blem is  somewhat  more  difficult*   It  involves  specifying  the  poles  and  zeros 
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in  such  a  manner  that  overshoot  can  be  expressed  in  terms  of  two  parameters 
instead  of  three.   Offhand  there  is  no  direct  procedure  but  rather  it  is 
one  of  trial  and  error.   The  following  examples  should  clarify  what  is  being 
attempted. 

1.    For  continuous  second  order  systems  the  complex  poles  can  be 


expressed  as  0"+  jw  or  w  (  £  +  j  y/  -  £2   ).  When  the  poles 
are  expressed  by  £  and  w  then  overshoot  becomes  a  function 

of  £  alone. 

2   2 
2.    Suppose  f(oCt3>^)-°^-/3      -  oL£+    Y  •   It  is  obvious 

that  f  can  be  expressed  as  a  function  of  two  variables  by  de- 
fining %  =  oC/3   .   Then 

f  =<<*,  p,  y  )  =  s2  -  s  +  r 

For  the  second  order  sampled  data  system  of  interest  a  successful  choice  of 
variables  has  been  found.   Define: 

Expression  of  c(nT)  in  terms  of  these  parameters  is  most  easily  accomplished 
by  recourse  to  geometry. 


Im[z] 


Re[z] 

Zl       1 
Figure   11 

Consider  figure  11.   By  inspection  the  folldwing  definitions  can  be  made, 
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^  =  argument  (p.-l) 

Vg=  argument  (P^^) 

<p  =  argument  (p.,) 

?-  arg  (Pl-Zl)  -  arg  (p^-1)  +  f  a  ^  -  ?,  4-  -J 
These  definitions  of   £  and  (^  are  consistent  with  those  made  earlier. 
An  interesting  representation  for  the  secant  of  £  results  from  these 
definitions.   Define  the  secant  in  terms  of  the  legs  of  right  triangle  (p  , 


zl$   b). 


sec  £  — 


ff-z. 


z.-b 


Now  observe  that  the  area  of  triangle  (p  ,  z  ,  1)  can  be  expressed  in  two 
ways 

area  A  (p  .  Zj.l)  -  {    (Uf>,-  ?,  \  \l  -'  Z,|)  =  -4  fl2'"  H  |l  -  R 


Vb 


Therefore 


2  \2 

=  _L  k-^l|/-z,| 


neyLg  Uz£lj£z5J  -  g  MlllzftJL      I5-Z« 
1  l^ll'-^l   '  li-zj    ■*|/?->||/-^| 


sec 


*-»*!?»-=* 


fc-')(P,-S) 


But  this  expression  is  identical  to  that  which  appears  in  the  equation  for 
c(nT).  Hence  c(nT)  can  be  written  as: 

c(nT)  =  |seC  £|  |P,  |  COs(/V$+f-7r) 
The  previously  defined  representation  of  "p."  as  a  function  of  £  and  w 
allows  expression  of  c(nT)  as: 

c(nT)  =|sec  £\  e  cos(//-£?  14  nT+£-ir) 

At  this  point  it  becomes  a  simple  matter  to  fit  a  continuous  curve  through 
the  discrete  output  c(nT).   To  do  this  simply  substitute  the  continuous 
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I 


variable  t  for  the  discrete  variable  nT. 

e.(t)  =  Isec  i\  e~SUJ"f  co5(i//-5£  W„+  +  S-t)  <59) 

c  (t)  is  the  continuous  function  which  will  be  used  in  place  of  the  true 
a 

output  c(t).   While  c  (t)  coincides  with  c(t)  at  sampling  instants  there 

a 

is  no  requirement  that  it  coincide  during  the  intersampling  periods  as 

well.   Nevertheless  it  is  assumed  that  c  (t)  is  "close  to"  c(t)  at  all 

times  and  that  the  system  output  is  adequately  represented  by  c  (t).  Under 

this  assumption  c  (t)  can  be  differentiated  to  obtain  peak  time,  maximum 

overshoot  and  other  characteristics  inherent  in  c(t). 

The  peak  time  is  found  by  differentiating  c  (t)  with  respect  to  t  and 

a 

equating  the  resulting  expression  to  zero. 

I 


T 
max 


^i^v      |f'-* 


Wt,,,  ^S|-S  +  7T 


(60) 


The  maximum  overshoot  is  found  by  substituting  the  expression  for  T    into 


c  (t). 

a 


c 
a 


-^[Mw-W-H, 


(max)  -.tf-'F  \Sec  £\  6      "FV!    L        U7r?I/      J!+   I  (61) 


The  portion  of  this  expression  underlined  is  the  maximum  overshoot  to  a 
unit  step  input.   Charts  of  overshoot  and  — ***     —   as  functions  of  £  and 
C      can  be  found  on  pages  175  and  176  of  Kuo  [2]. 

The  equations  for  T    (60)  and  c  (max),  (61)  characterize  a  second 

Hlcl«v  ex 

order  sampled  data  system  in  terms  of  two  parameters.   However  they  are  not 
exact  representations  and  the  problem  is  not  really  complete  until  regions 
of  validity  for  this  characterization  are  determined.   Lindorff  (by  what 
seems  to  be  intuitive  argument)  has  determined  that  the  approximation  is 
good  if  the  complex  poles  are  in  the  first  and  fourth  quadrants  and  are  not 
too  close  to  the  origin  of  the  z-plane  [9].   [Poles  in  the  first  and  fourth 
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quadrants  insure  at  least  four  samples  per  cycle  of  the  dominant  oscilla- 
tory mode  and  thus  give  a  reasonably  good  approximation  for  this  mode.  How- 
ever, if  the  poles  are  too  close  to  the  origin  the  rapid  decay  offered  by 
the  Ipl   factor  will  predominate  to  the  point  where  the  z-transform  [and  thus 

the  approximating  function  c  (t)]  no  longer  adequately  characterize  the 

a 

system.   That  is,  the  sampling  rate  is  too  slow  to  follow  the  decay]. 

Analysis  of  Dominance  Requirements 

For  the  z-plane  dominant  mode  model  developed  in  Kuo  and  reviewed  in  the 
preceding  section  to  have  value  the  conditions  under  which  model-like  be- 
havior prevails  must  be  determined.  These  requirements  can  be  established 
by  considering  the  inverse  z-transform  integral. 


*1J 


c(nT)  =  -L-   I    CCDl"-'    dl 


r 
where  P  is  a  path  which  encloses  all  singularities  of  C(z).  But 

C(z)  =  G(z)R(z) 

C(z)  ■  ]J  CZ-Zk)  — Z__  for  a  unit  step  input. 

jf(z-r,)  z-i 

Substituting  this  expression  into  the  integral  equation  yields 

i    r  it  (z-z»)     z   ,*j, 

c(nT)  -  — — -  /  — *r 2  dZ 

2irW  ft(z-5)     z-i 

This  integral  can  now  be  evaluated  by  the  theory  of  residues.   Expressed  in 
this  manner  it  becomes: 

c(nT)  -  >   5=   1!  ^~Zk)   ,  U  V~^k)  ^r steady  state  output 

£*-*£(&- 6)       n(i-pj" 

Analysis  of  this  expression  reveals  the  two  conditions  necessary  for  dominant 
mode  behavior. 

1.   All  but  the  two  dominant  poles  must  lie  near  the  origin  of 
the  z-plane.   After  a  small  number  of  samples  the  residues 
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contributed  by  poles  near  the  origin  are  negligible;  the  pn 
term  dominates  these  residues  and  as  a  factor  it  becomes 
almost  zero. 
2.   Excess  zeros  over  the  real  zero  employed  in  the  second  order 
model  must  lie  near  the  origin  of  the  z-plane.   When  so 
located  two  effects  are  assumed  to  result. 

a.  The  residues  of  poles  near  the  origin  (which  are 
neglected  in  the  second  order  model)  will  be  small. 

b.  These  zeros  form  dipoles  (pole-zero  pairs  close  to- 
gether) which  have  almost  negligible  effect  on  the 
residue  of  the  dominant  poles,  [i.e.   If  (p . -*.  )  — 
(p  -p,  )  then  these  common  factors  cancel  in  the  ex- 
pression for  the  dominant  pole  p . ] .   (Where  time  de- 
lays are  absent,  dominant  mode  design  requires  as 
many  poles  as  zeros  to  be  ignored  or  lie  near  the 
origin.   This  proximity  of  ignored  poles  and  ignored 
zeros  essentially  guarantees  dipole  formation  when 
viewed  from  the  location  of  the  dominant  poles) . 

Principles  of  Dominant  Mode  Compensation  in  the  Parameter  Plane 

With  model  behavior  and  restrictions  to  insure  model  behavior  in  higher 
order  system  established  attention  turns  to  the  problem  of  properly  compensat- 
ing existing  systems  to  resemble  the  model.   Parameter  plane  methods  are  well 
suited  to  this  aspect  of  dominant  modeldesign.   This  approach  not  only  allows 
selection  of  the  dominant  (model)  poles  but  also  guarantees  dominant  mode  (or 
model  like)  behavior  by  forcing  compliance  with  the  restrictions  mentioned 
above.  Tacit  is  the  assumption  that  parameter  values  (  oC   ,  /3  )  exist  which 
allow  accomplishment  of  this  ambitious  objective.   Justification  of  these 

78 


statements  can  be  had  by  recourse  to  Chapter  III.   Recall  that  the  image 
of  a  z-plane  contour  divides  the  parameter  plane  into  regions  whose  distin- 
guishing characteristic  is  that  any  parameter  setting  within  a  given  region 
places  an  equal  number  of  characteristic  equation  roots  within  the  z-plane 
contour  which  was  mapped.   Suppose  a  circular  contour  centered  at  the  origin 
of  the  z-plane  and  having  a  radius  small  enough  so  that  all  poles  within 
it  can  be  considered  negligible,  is  mapped  through  the  characteristic  equa- 
tion into  the  parameter  plane.   Of  the  parameter  plane  regions  defined  by 
this  mapping  the  dominant  mode  designer  considers  only  those  for  which  all 
roots  but  two  of  the  characteristic  equation  lie  within  the  z-plane  contour. 
If  within  such  a  region  parameter  settings  exist  which  place  two  poles  near 
the  unit  circle  then  these  poles  must  be  dominant  (because  it  is  known  that 
no  other  poles  except  those  near  the  origin  exist) .   To  find  these  poles 

auxiliary  contours  of  constant  w  and  constant  C     are  mapped  through  the 

n  ^ 

characteristic  equation  into  the  parameter  plane.  The  intersection  of  an 
acceptable  w  contour  and  t>    contour  in  the  region  mentioned  above  constitu- 
tes a  set  of  parameter  values  which  simultaneously  guarantees  dominance  and 
the  desired  pole  locations. 

To  this  point  the  closed  loop  zero  specification  required  by  the  domin- 
ant mode  model  has  been  ignored.   Such  slight  is  intentional  since  zeros  of 
the  open  and  closed  loop  transfer  function  are  identical  and  are  not  general- 
ly subject  to  modification  by  available  system  parameters.   Therefore  dis- 
counting cancellation,  zeros  should  be  considered  as  entities  which  allow  or 
disallow  dominant  mode  compensation  rather  than  as  entities  to  be  modified. 

Dominant  mode  design  using  parameter  plane  concepts  is  readily  extended 
to  continuous  systems  described  in  the  s-plane.   Instead  of  a  small  circle 
near  the  origin  a  contour  of  constant  0"  (  6~  =C  w  )  sufficiently  removed  from 
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the  imaginary  axis  is  mapped  through  the  characteristic  equation.   The  image 
of  this  s-plane  contour  divides  the  parameter  plane  in  exactly  the  same 
manner  as  the  image  of  the  small  circle  in  the  z-plane;  subsequent  s  and 
z-plane  procedures  are  identical.   In  fact  dominant  mode  design  in  the  s- 
plane  is  generally  simpler  than  in  the  z-plane  because  the  s-plane  dominant 
mode  model  has  no  zeros  in  the  numerator  of  its  transfer  function. 

Cascade  Compensation  (for  Dominant  Mode  Design) 

The  principles  discussed  in  the  preceding  chapter  fairly  well  dictate 
the  form  that  a  cascade  compensator  must  have  if  dominant  mode  design  is 
to  be  achieved.   Consider  the  sampled  data  system  of  figure  (12). 


R(z) 


Figure  12 
G(z)  is  the  plant  transfer  function  which  is  assumed  to  be  of  the  form: 

«.)■  k/!(Z  +  ^  (62) 

/I  (z+£) 

The  system  open  loop  zeros  provide  the  best  clue  as  to  the  feasibility  of 
dominant  mode  design.   Since  the  zeros  of  the  open  and  closed  loop  transfer 
functions  are  identical,  all  zeros  except  one  must  lie  near  the  origin  of 
the  z-plane.  When  this  condition  is  not  met  no  cascade  compensation  is 
capable  of  driving  the  system  to  the  second  order  dominant  mode  desired. 
It  is  also  prudent  to  observe  the  system  poles  before  undertaking  dominant 
mode  design.   Cascade  compensation  can  succeed,  disallowing  pole-zero  can- 
cellation, only  if  addition  of  poles  and  zeros  offers  the  possibility  of 
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reshaping  the  root  locus  so  that  all  but  two  complex  poles  and  one  real  zero 
lie  near  the  origin.   If  the  general  uncompensated  system  pole-zero  configura- 
tion is  such  that  this  is  manifestly  impossible  then  there  is  no  point  in 
employing  parameter  plane  procedures  to  prove  it. 

These  observations  can  be  made  more  rigorous  by  translating  them  into 
restriction  on  the  transfer  function  of  the  cascade  compensator. 

1.  The  zeros  of  the  cascade  compensator  must  lie  near  the  origin 
of  the  z-plane  (assuming  that  one  system  zero  is  not  so 
located).   In  this  way  these  zeros  can  be  neglected  without 
affecting  the  validity  of  the  dominant  second  order  approxi- 
mation. 

2.  The  numerator  should  be  of  the  same  order  in  z  as  the  de- 
nominator (primarily  because  time  delay  is  an  undesirable 
property  in  a  compensator) . 

3.  The  steady  state  gain  of  the  compensator  should  be  unity. 
With  these  properties  in  mind  the  transfer  function  of  the  cascade  compen- 
sator must  have  the  form: 


Q  +  a^ri 4-Q,  -h  Cl0)7:n 


•.<•>  -  V+a,Mz~+.:.+al2  +  a,  (63) 

With  this  cascade  compensator  installed  the  closed  loop  transfer  function 
of  the  system  in  figure  (12)  is  expressible  as: 


M 

C(z)  . k  (n-o*-,  +  •  -  ■  +  a,+  a.)  z"  l\  (Z~ZQ 


R(z) 


(z'+a^z"-'  +■   -+ao)/|  (z-eL)+k(^^,+  -+a0)z"  /Kz-z.) 


(64) 


t 


where  k,  z  ,  and  p.  are  known  fixed  parameters  off the  plant.   The  a.'s  are 
now  chosen  to  yield  a  satisfactory  dominant  mode  second  order  system.  While 
equations  (63)  and  (64)  appear  formidable,  parameter  plane  methods  will  yield 
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a  solution  almost  by  Inspection  especially  if  but  two  of  the  a.'s  are 
variable. 

This  method  is  also  applicable  when  pole-zero  cancellation  is  allowed, 
In  this  case  the  cascade  compensator  has  a  transfer  function  of  the  form: 

At 

G  (z)  =  '(!+(***+■■   +-a,  t- cpz"  il,(z"^ 

where  (z  +  r.)  and  (z  +  t  )  cancel  poles  and  zeros  of  the  plant  transfer 
function  respectively. 


R(z) 


Feedback  Compensation  (for  Dominant  Mode  Design) 

Feedback  compensation  instead  of  cascade  compensation  can  be  used  to 
effect  dominant  mode  design.   As  with  cascade  compensation  certain  restric- 
tions are  placed  on  the  feedback  transfer  function  if  the  chosen  second 
order  dominant  model  is  to  be  valid.   Consider  the  sampled  data  system  of 
figure  (13). 

where  H  =  numerator  of  H(z) 

H  =  denominator  of  H(z) 
G  ■  numerator  of  G(z) 
G  ■  denominator  of  G(Z) 


Figure  (13) 
Assume  that: 

a.  The  plant  transfer  function  G(s)  contains  at  least  one  integra- 

tion,  (i.e.   The  denominator  of  G(s)  contains  a  factor  of  s 
or  the  denominator  of  G(z)  contains  a  factor  [z-lf). 

b.  At  most  all  but  one  zero  of  the  plant  transfer  function  is 

located  near  the  origin  of  the  z-plane. 
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Requirements: 

1.  Design  a, feedback  compensator  H(z)  which  yields  a  second  order 
dominant  mode  system  with  acceptable  time  response  characteristics. 

2.  The  system  must  follow  a  step  input  with  zero  steady  state  error. 
The  first  step  toward  solution  will  be  to  apply  the  assumptions  with  the 
objective  of  obtaining  an  acceptable  class  of  transfer  functions  H(z). 
Parameter  plane  analysis  applied  to  this  class  will  then  yield  acceptable 
H(z)'s  if  they  exist.   The  closed  loop  transfer  function  for  the  system  of 
figure  (13)  is 

R(z)        /'+  G(z)H(z) 

C(Z)      G>"Ud 
R(z)        GdHd+G"H" 

The  first  restriction  on  H(z)  can  be  obtained  by  observing  that  poles  of 
H(z)  are  zeros  of  the  closed  loop  transfer  function.   To  effect  dominant 
mode  design  these  closed  loop  zeros  [or  equivalently  poles  of  H(z) ]  must  lie 
near  the  origin.   The  second  restriction  on  H(z)  results  from  the  requirement 
that  there  be  zero  steady  state  error  to  a  step  input.   This  restriction  is 
obtained  by  applying  the  final  value  theorem  to  the  closed  loop  transfer 
function. 


for  unit  step 
input 


(       '     z-i  [_  z  (z^ij  GJH7+&"H"J   U"H  +G"H  J 

But  z-w  [G  ]  ■  0  since  the  plant  transfer  function  was  assumed  to  contain  an 
integration.   Hence: 
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or 


lim  __r  x    . 
z^/  H(z)  -  1 


To  meet  these  requirements  the  feedback  path  must  have  a  transfer  function 
of  the  form  H 

(z-/)J  aKzK 

H(z)  -  1  + &E& (65) 

This  transfer  function  is  slightly  more  restrictive  than  necessary  since  all 
poles  of  H(z)  were  placed  at  the  origin  (for  computational  ease)  rather 
than  near  the  origin  which  is  the  specific  requirement.  Also  notice  that 
H(z)  is  written  as  containing  no  time  delays  although  this  is  not  a  specific 
requirement.  Time  delays  can  easily  be  added  by  increasing  the  powers  of 
z  in  the  denominator  of  H(z) . 

Example  of  Dominant  Mode  Design  in  the  Parameter  Plane 

A  numerical  example  is  appropriate  to  illustrate  theory  and  to  rein- 
force stated  results.   Suppose  the  plant  shown  in  figure  14  is  to  be  compen- 
sated to  a  dominant  mode  which  satisfies  the  specifications 

i,2<IJ  <  1.4 
^   pt  ^ 

rise  time   -   less  than  5  sampling  periods 
As  an  initial  step  in  the  compensation  procedures  the  zeros  of  this  plant's 
closed  loop  transfer  function  (equation  66)  are  checked  for  compliance  with 
necessary  conditions  for  dominant  mode  design.   Namely,  validity  of  dominant 
mode  design  presupposes  all  zeros  but  one  are  sufficiently  near  the  origin 
to  be  negligible.   The  real  zero  at  z  «  -.295  stretches  this  requirement  but 
when  compared  to  the  model  zero  at  z  =  -3.72  it  will  probably  not  be  suf- 
ficient to  destroy  dominant  mode  behavior.   Tables  describing  model  be- 
havior such  as  those  referenced  in  the  first  section  of  this  chapter  are  now 
consulted  for  £,  £  values  which  satisfy  the  specifications.  For  this 
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The  closed  loop  transfer  function  of  the  system  shown  above 


is: 


C(z)   N(z) 


R(z)   D(z) 


where 


N(z)  =  0.022(z2)(z+3.72)(z+0.29£)(z+0.0l£) 

D(z)  =  z6  +  (0.0220/S  ^1.8llj.)z£ 

+  (0.0220* +O.0667/?  +1.066)  z^ 
+  (0.o667«<-o.o632/« -0.1l6)z3 
+  (-0.0632o<  -0.0252/S  4-0.00128)z2 
4-  (-0.02£2<< -0.00037/*  )z 
+  (-0.00037 °<  ) 


(66) 
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example  approximate  value  ranges 

-30°  <  S  <  10° 

0.3  ^  £  <   0.45 
appear  satisfactory. 

The  parameter  plane  is  now  employed  to  effect  dominant  mode  behavior  in 

the  actual  system.   Those  system  poles  whose  settling  time  is  five  times  or 

more  faster  than  that  of  the  dominant  mode  will  be  considered  negligible. 

Since  w    (the  settling  time  indicator)  of  the  dominant  mode  is  about  unity 
nz 

all  other  poles  must  lie  within  a  circle  of  radius  w    =0.2  centered  at 

«z 

the  origin.   The  image  of  this  settling  time  contour  divides  the  parameter 
plane  into  regions  containing  an  equal  number  of  roots  to  the  characteris- 
tic equation  (figure  15);  the  region  for  which  all  but  two  characteristic 
roots  lie  within  the  designated  settling  time  contour  is  shaded.  The  pre- 
sence of  constant  C,     curves  in  this  region  within  the  range  needed  to  meet 
specifications  is  ascertained  and  on  each  curve  the  range  of  satisfactory  d 
values  encountered  in  traversing  the  region  is  marked.   These  actions  serve 
to  delineate  a  parameter  plane  region  (shaded  in  figure  16)  within  which  all 
parameter  settings  provide  a  second  order  dominant  system  meeting  specifica- 
tions.  Figure  17  compares  the  transient  response  of  a  true  second  order 
system  with  the  response  of  the  system  in  figure  14  compensated  for  dominant 
mode  behavior.  The  dominant  poles  and  zeros  in  the  latter  case  identify 
with  the  poles  and  zero  of  the  true  second  order  system  whose  response  is 
displayed.   Numerical  values  for  the  compensated  system  are  as  follows: 

parameter  settings 

oC  m   -0.13  /S  -  0.73 

dominant  closed  loop  poles  and  zero 

poles:  z  =  0.776  +  jO.362         zero:  =  z  ■  -3.72 
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ignored  closed  loop  poles  and  zeros 

poles: i  z  »  0.182  +  jO.078        zeros:  z  ■  -0.295 

z  -  -0.016  z  -  -0.015 

2 
z  -  -0.106  z  ■  0.0 
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CHAPTER  VI 

Frequently  the  niceties  of  a  compensation  procedure  are  more  clearly  revealed 
by  numerical  examples.   In  deference  to  this  hypothesis  several  examples 
emphasizing  various  aspects  of  parameter  plane  theory  are  presented  in  this 
chapter. 

Example  of  Compensation  by  Parameter  Plane  Technique 

The  fragmentary  examples  in  preceding  chapters  relating  to  the  system 
shown  in  figure  (2)  can  be  gathered  into  a  fairly  detailed  and  representative 
illustration  of  parameter  plane  technique.   Typically  the  engineer  is  given 
a  basic  plant  such  as  that  shown  in  figure  (2)  and  he  is  told  to  modify  its 
performance  to  meet  specification.   Suppose  these  specifications  require 
maximum  overshoot  to  step  inputs  to  be  under  twenty  per  cent,  limit  the 
sampling  rate  to  two  samples  per  second,  and  require  steady  state  errors 
to  unit  ramp  inputs  of  less  than  one.   In  addition  short  rise  time  and  short 
settling  time  is  desirable.  Direct  consideration  of  all  these  specifications 
is  not  practical  in  the  parameter  plane  but  tractable  indicators  for  each 
exist.   These  specifications  modified  for  parameter  plane  utilization  be- 
come: 

!•   %-«»«>  °-5  l»UeuofMp(:  <l-2 

2.  E     <T  1   direct  consideration 

s.s.  ^ 

3.  large  bandwidth    in  lieu  of  short  rise  time 

4.  short  settlint  time   direct  consideration 

The  engineer  decides  that  if  feasible  a  combination  of  gain  adjustment  and 
derivative  feedback  offers  the  most  attractive  way  to  meet  these  specifica- 
tions.  The  closed  loop  z- transfer  function  with  the  perspective  compensa- 
tion incorporated  (equation  32)  is  amenable  to  parameter  plane  analysis. 
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Pertinent  curves  Sufficient  to  determine  the  parameter  settings  required  to 
effect  a  satisfactory  system  are  shown  in  figure  18.   One  observes  that  to 
meet  overshoot  specifications  parameter  settings  ( oL  ,  /6  )  must  lie  to  the 
left  of  the  £  ■  0.5  contour  while  steady  state  error  requires  these  set- 
tings to  be  to  the  right  of  the  error  locus  shown.  Unfortunately  the  inter- 
section of  these  two  areas  is  rather  small  leaving  the  designer  little  room 
to  maneuver  for  the  desirable  properties  of  short  settling  time  and  short 
rise  time.   Nevertheless  a  constant  settling  time  curve  is  shown;  all  para- 
meter settings  within  the  contour  yield  a  system  with  a  shorter  settling 
time  than  is  obtained  for  parameter  settings  outside  the  contour.   Band- 
width curves  which  are  indicative  of  rise  time  are  also  drawn;  regions  to 
the  right  of  each  curve  have  shorter  rise  times  than  those  to  the  left. 
From  these  parameter  plane  considerations  one  concludes  that  settings  of 
o4=  k   -  0.5 
/$»  kkt  ■  0.5 
should  yield  an  acceptable  system.   The  time  response  (at  the  sampling 
instants)  for  this  choice  of  parameters  is  displayed  for  step  and  ramp  in- 
puts in  figure  19a  and  19b  respectively. 

Example  of  Analysis  by  Parameter  Plane  Technique 

In  addition  to  its  capabilities  as  a  design  tool  the  parameter  plane 
also  provides  an  excellent  vehicle  for  analysis.  Consider  the  second  and 
third  order  feedback  systems  shown  in  figure  20a  and  20  b  respectively  and 
suppose  them  to  be  components  of  a  still  larger  system.   If  pole  b  is  far 
removed  from  pole  a  along  the  negative  real  axis  it  is  well  known  that  the 
third  order  system  will  behave  essentially  like  its  second  order  counterpart. 
However  approximate  second  order  behavior  of  this  component  depends  not  only 
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on  the  distance  that  pole  b  is  "removed"  but  additionally  it  depends  on 
gain  values  k  and  kfc.  A  knowledge  of  when  second  order  behavior  prevails 
as  a  function  of  k,  kfc,  and  b  is  of  consequence  in  analyzing  the  overall 
system.  Furthermore,  if  k  and  kft  are  used  to  adjust  overall  system  perfor- 
mance -  for  what  range  of  values  can  these  adjustments  be  treated  as  ad- 
justments to  a  second  order  system  and  when  must  the  third  pole  be  consider- 
ed? These  questions  are  conveniently  answered  on  the  parameter  plane, 
shown  in  figure  21,  where  it  is  assumed  that  nominally  these  poles  a#e  ad- 
justed to  a  zeta  value  of  0.5.   (If  additional  £  values  are  of  interest 
these  contours  are  also  mapped  into  the  parameter  plane).   Proximity  of 
second  and  third  order   £  =  0.5  curves  denotes  approximately  equal  charac- 
teristic equation  roots  for  parameter  settings  in  these  regions;  consequent- 
ly the  third  order  system  is  closely  approximated  by  its  second  order  count- 
erpart. Figure  21  shows  such  approximations  to  be  valid  for 

0.6  <o(=k   <  1.6 

for  b  =  10 

0  <(S  ±   kkt<   0.6 

0.6  <  OL      <  1.2 

for  b  =  5 

0  <   /3   <0.35 
Notice  that  even  a  10:1  ratio  between  poles  a  and  b  does  not  justify  second 
order  approximation  f or  *>  2.4  or  p  >  1.5.  In  figure  22  the  effect  of 
these  approximations  on  bandwidth  is  shown. 

The  author  cautions  against  haphazard  extension  of  these  approxima- 
tions. For  instance,  proximity  of  £  -  curves  for  a  second  and  fourth 
order  system  would  not  be  sufficient  to  guarantee  the  validity  of  a  second 
order  approximation.  Not  only  must  the  location  of  the  other  two  roots  to 
the  fourth  order  system  be  investigated  but  the  region  of  proximity  for  the 
£  -  curves  must  be  checked  for  a  corresponding  proximity  of  wn  values  as 

well. 
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Advanced  Example  of  Compensation  by  Parameter  Plane  Technique 

As  mentioned  in  the  introduction  most  graphical  methods  are  seriously 
limited  when  called  upon  to  display  system  performance  measures  as  a  func- 
tion of  more  than  one  variable.   In  this  respect  parameter  plane  methods 
have  definite  advantages  for  performance  indicators  such  as  damping,  settling 
time  and  bandwidth  are  inherently  displayed  as  functions  of  two  system  vari- 
ables.  Furthermore,  when  only  a  single  zeta  value  (or  a  single  settling 
time  or  a  single  bandwidth)  is  of  interest  parameter  plane  methods  compet- 
ently display  the  effect  of  three  variables  on  system  performance.   Even 
more  variables  can  be  considered  without  undue  effort  as  will  be  illustrat- 
ed in  designing  compensation  for  the  system  of  figure  23.  Here  the  design- 
er has  four  parameters  at  his  disposal,  the  system  gain,  the  amount  of  de- 
rivative feedback,  and  the  pole  and  zero  locations  of  the  single  section 
filter.   Settings  for  these  parameters  which  allow  the  system  to  meet  the 
following  performance  specifications  are  to  be  determined. 

1.  Maximum  overshoot  must  be  less  than  1.25  or  zeta  is  to  be 
equal  to  or  greater  than  0.4. 

2.  The  steady  state  error  to  unit  ramp  inputs  must  be  less 
than  0.25. 

3.  The  bandwidth  (-3  db  point)  of  this  system  must  occur  at  a 
radian  frequency  greater  than  7T  to  allow  sufficiently  fast 
response. 

Examination  of  the  closed  loop  transfer  function  reveals  its  suitability 
for  parameter  plane  analysis  with  oC  and  /3    capable  of  definition  in  any 
of  several  ways.  The  curves  illustrated  result  from  the  choice 

oC=k,      /?=kp/q 
with  p  and  k  considered  as  third  and  fourth  parameters.   Initially  k  is 
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equated  to  zero  and  the  efficacy  of  lag  and  lead  compensation  is  investi- 
gated.  The  effect  of  such  compensation  on  stability  and  bandwidth,  as 
displayed  in  a  parameter  plane  format,  is  shown  in  figures  24  and  25  res- 
pectively.  That  neither  lag  nor  lead  compensation  is  capable  of  rendering 
a  satisfactory  system  is  immediately  evident.   The  intersection  of  the  para- 
meter plane  region  denoting  stability  with  that  region  providing  a  satis- 
factory steady  state  error  is  empty  for  lead  compensation  whereas  the  inter- 
section of  those  regions  denoting  stability  with  those  denoting  acceptable 
bandwidth  is  empty  for  lag  compensation.  However  the  parameter  plane  dis- 
play does  indicate  that  a  "good  location"  for  a  lead  section  pole  is  around 
s  ■  -8  and  that  a  lag  pole  might  profitably  be  located  at  s  ■  -0.05. 

The  effect  of  derivative  feedback  in  combination  with  lead  and  lag 
filters  can  now  be  investigated.  The  damping  curves  of  figure  26  and  the 
bandwidth  curves  of  figure  27  display  the  results  for  various  values  of  de- 
rivative feedback  used  in  conjunction  with  a  lead  filter  (or  more  accurately 
a  single  section  filter  whose  pole  is  located  at  s  ■  -8).   Figure  28  dis- 
plays an  analogous  result  for  feedback  used  in  conjunction  with  a  single 
section  filter  whose  pole  is  located  at  s  =  -0.05. 

In  both  instances  system  performance  is  enhanced  by  employment  of  some 
derivative  feedback.   In  this  case  lag  filters  provide  the  easiest  way  to 
improve  steady  state  accuracy  but  even  with  feedback  their  bandwidth  limita- 
tions (slow  rise  times)  cannot  be  overcome.   On  the  other  hand  lead  filters 
readily  provide  acceptable  bandwidth  and  with  derivative  feedback  both  damp- 
ing and  error  specifications  can  be  met.   Intersection  of  acceptable  damp- 
ing, error,  and  bandwidth  regions  indicate  parameter  values  of 

k  =  oL    =  5.0  kfc  =  0.5 


should  suffice. 


p  ■  8.0  q  =   -^-p  «  3.0 
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All  curves  shown  are  images  of  the  unit 
circle  in  the  z-plane  and  represent  the 
stability  limit  for  parameter  settings 
«C  and/Q.  For  each  "p""  shown  the  stable 
region  is  to  the  left  of  the  image  curve 
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FIGURE  25 

THE  EFFECT  OF  THE  POLE  LOCATION  OF  A  SINGLE  SECTION 
COMPENSATOR  ON  SYSTEM  BANDWIDTH  AND  STEADY  STATE  ERROR 
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Curves  show  the  image  of  a  constant  damping  spiral 
in  the  z-plane,  equivalent  to  £=0*\\.9    as  affected 
by  varying  amounts  of  derivative  feedback. 
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EFFECT  OF  DERIVATIVE  FEEDBACK 
ON  A  LEAD  COMPENSATED  SYSTEM 
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All  bandwidth  loci  are  specified  as 
A—  -3  db   at  frequency  w=TT 
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and  the  various  curves  represent 
2  the  effect  of  varying  the  amount  of 
|  derivative  feedback.   Regions  above, 
{  to  the  right  or  located  within  a 
j  bandwidth  locus  have  bandwidths 
'  greater  than  that  indicated  by 
parameter  settings  on  the  locus. 
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CHAPTER  VII 

The  general  availability  of  the  digital  computer  has  fostered  revolu- 
tionary advances  in  servo-system  and  especially,  in  sampled  data  system 
design.   However  many  of  these  new  approaches  completely  discard  rather 
than  refurbish  existing  methods  and  thus  void  much  of  the  experience  gained 
by  practising  engineers.   In  contrast  parameter  plane  methods  retain  the 
principle  of  design  by  such  familiar  specifications  as  damping  and  settl- 
ing time  but  where  classical  methods  could  display  the  effect  on  only  one 
variable  on  these  specifications,  parameter  plane  methods  can  display  the 
simultaneous  effect  of  two.  This  paper  has  extended  parameter  plane  methods 
in  two  directions.   First  the  requirement  that  parameters  appear  linearly 
in  the  characteristic  equation  was  eased  to  permit  consideration  of  those 
cases  in  which  their  product  is  also  present.   Extension  in  the  second 
direction  increased  the  information  which  could  be  displayed,  and  hence  de- 
signed for,  on  the  parameter  plane  to  include  bandwidth  and  steady  state 
error  specifications.   This  latter  extension  is  particularly  significant 
because  it  relates  performance  to  the  zeros  of  the  closed  loop  transfer 
function  which  were  previously  ignored.  Another  novel  capability  of  the 
parameter  plane  which  was  introduced  in  Chapter  5  is  the  simultaneous  de- 
sign of  a  dominant  mode  system  and  the  guarantee  of  such  dominance. 

A  great  deal  of  further  investigation  into  parameter  plane  methods  in 
both  theory  and  application  still  remains.   To  amplify  questions  typical  of 
those  to  be  resolved  consider  the  consequences  of  allowing  product  terms  to 
appear  in  the  coefficients.  Now  points  in  the  s-  or  z-plane  no  longer  have 
unique  images  in  the  parameter  plane  and  some  points  in  the  transform  planes 
could  conceivably  have  no  real  image  at  all.   It  is  easy  enough  to  postulate 
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that  available  parameters  cannot  place  characteristic  equation  roots  at 
these  latter  locations  and  that  closed  contours  in  the  transform  planes  will 
result  in  consistent  contours  in  the  parameter  plane.   That  is  -  these  para- 
meter plane  contours  are  expected  to  be  consistent  with  the  interpretation 
expounded  in  Chapter  3.  While  such  postulation  seems  reasonable  enough 
for  the  product  case,  it  certainly  warrants  a  more  rigorous  framework  if 
these  methods  are  to  be  extended  still  further. 

To  the  author's  knowledge  the  bandwidth  and  error  curves  represent  the 
first  display  of  information  on  the  parameter  plane  which  is  not  directly 
obtainable  from  the  characteristic  equation.   Research  to  discover  other 
performance  indicators  capable  of  being  displayed  on  the  parameter  plane 
should  certainly  be  fruitful.   In  fact  the  author  has  obtained  the  equa- 
tion governing  loci  of  constant  root  sensitivity,  but  each  point  thereon 
must  be  computed  by  an  iterative  procedure  and  it  is  too  unwieldy  to  be 
generally  useful.  Utilization  of  this  auxiliary  information  in  conjunction 
with  that  normally  available  on  the  parameter  plane  is  illustrated  in 
Chapter  6.  Another  facet  of  parameter  plane  methods  illustrated  by  these 
problems  is  the  relative  ease  with  which  the  effect  of  a  third  variable 
parameter  can  be  analysed.   This  ease  represents  a  significant  improvement 
over  the  almost  unworkable  situation  a  three  variable  problem  presents  to 
root  locus  methods. 

The  author  feels  that  parameter  plane  methods  as  enhanced  by  this  paper 
should  equal  or  exceed,  for  two  parameter  problems,  those  capabilities  which 
the  root  locus  brings  to  problems  of  one  parameter.   Furthermore,  the  re- 
liance on  familiar  performance  indicators  minimizes  prerequisite  studies 
and  enhances  the  experience  of  most  engineers.   These  traitJs  auger  its  more 
widespread  appearance. 
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APPENDIX 

Some  Self -Explanatory  Parameter  Plane  Curves 

These  self-explanatory  curves  should  serve  to  relate  the  familiar 
second  order  system  to  some  parameter  plane  representations  of  it. 
A  knowledge  of  the  general  pattern  of  such  relations  will  serve  the 
reader  when  interpreting  parameter  plane  curves  for  more  complex 
systems. 
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